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CHAPTER  I 


INTRODUCTION 

The  transmission  of  sound  in  the  ocean  has  been  studied  for  many 
years.  These  studies  have  frequently  been  directed  toward  the 
measurement  and  analysis  of  the  ways  in  which  the  sound  waves  are 
attenuated  as  they  propagate  out  from  the  source  of  acoustic  energy. 

In  spite  of  this  effort,  many  problems  remain  unsolved.  This  study 
attempts  to  expand  the  body  of  knowledge  related  to  one  aspect  of  the 
attenuation  problem. 

Normal-mode  theory  has  been  most  successful  in  describing 
qualitative  propagation  effects  for  range-invariant,  shallow-water 
environments.  Given  Che  relatively  simple  assumptions  of  the  theory, 
comparisons  of  experimental  and  calculated  propagation  loss  do  not 
always  correlate  well.  Weston  (1967)  recognized  this  problem  and 
suggested  such  reasons  as  boundary-related  losses  and  losses  of 
biological  origin.  Boundary  losses.  In  particular  those  related  to 
the  seabed,  are  the  primary  concern  in  this  study. 

The  work  of  Pekeris  (1948)  was  probably  the  first  attempt  to 
explain  shallow-water  propagation  in  terms  of  normal-mode  theory. 
Perturbations  of  the  basic  theory,  to  account  for  boundary  losses, 
were  introduced  by  Kornhauser  and  Raney  (1955),  and  Tolstoy  (1958). 

Most  recently,  the  theory  has  been  successfully  compared  with 
experimental  propagation  results  by  Kibblewhite  and  Denham  (1966), 
Ferris  (1972),  and  Ingenito  (1973). 


At-sea  measurements  and  model  tank  experiments  by  Eby  et  al. 
(1960)  have  demonstrated  the  validity  of  the  theory  and  its  modification 
for  lossy  bottoms. 

As  can  be  seen  from  this  discussion,  the  theory  and  its 
modification  for  lossy  bottoms  have  been  shown  to  work.  However, 
suppose  it  is  desired  to  predict  propagation  loss  for  an  area  where  a 
limited  quantity  of  bottom  information  such  as  bottom  type  (i.e.,  sand, 
mud,  etc.)  and  historical  sound-velocity  information  is  available. 

The  question  arises:  Is  it  possible  to  predict  modal-attenuation 
coefficients  and  sufficient  sub-bottom  Information  so  that  transmission 
loss  may  be  predicted  through  the  use  of  a  normal-mode  model?  The 
answer  to  this  question  is  heavily  dependent  upon  our  ability  to  predict 
bottom  absorption,  a  topic  to  be  addressed  next. 

Volume  absorption  in  the  sediment  is  generally  the  dominant  loss 
mechanism  for  flat  gravel,  sand,  or  mud  bottoms.  The  physical  basis  for 
this  loss  and  mathematical  models  describing  propagation  in  a  fluid- 
filled  porous  sediment  have  been  proposed  by  Morse  (1952),  Biot  (1956, 
1962),  Brutsaert  (1964),  and  Walton  (1977).  Generally,  two  sources  of 
dissipation  have  been  considered;  that  is,  viscous  losses  due  to  fluid 
movement  within  the  skeletal  frame  and  friction  losses  due  to  movement 
between  particles  within  the  frame  itself.  The  work  of  Biot  and,  more 
recently  following  his  theory,  the  work  of  Stoll  and  Bryan  (1970),  and 
Stoll  (1974,  1977,  1979)  have  found  the  greatest  acceptance  of  the 
mathematical  theories.  The  work  of  Hamilton  (1972,  1974),  mostly  of 
an  empirical  nature,  has  also  received  some  acceptance  for  predicting 
absorption  values. 
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In  this  study,  the  determination  of  volume  absorption  in  the 
sediment  is  approached  from  two  aspects.  First,  bottom  absorption  is 
obtained  as  a  function  of  frequency  from  transmission  loss  data  and 
calculated  mode-attenuation  coefficients.  Second,  using  data  from  grab- 
sample  and  core  analyses  as  input  to  the  sediment  models  of  Biot  and 
Hamilton,  bottom  absorption  is  calculated  as  a  function  of  frequency. 
After  establishing,  by  comparison  with  acoustic  data,  that  the  sediment 
model  can  be  used  to  predict  bottom  absorption,  the  sediment  model  and 
other  seabed  information  is  used  in  a  normal-mode  model  to  predict 
transmission  loss  in  test  areas  where  comparative  data  is  available. 

Without  adequate  experimental  data,  the  study  would  be  impossible. 
Fortunately,  propagation  data  is  available  for  experiments  performed 
off  the  Florida  and  Texas  coasts,  and  one  experiment  performed  on  the 
Scotian  Shelf.  These  sites,  in  general,  are  well-documented  with  regard 
to  the  sediments  and  sub-bottom  structure  present. 

Following  the  introductory  material  presented  in  this  chapter, 
a  discussion  of  normal-mode  theory  and  the  theory  relating  to 
propagation  in  fluid-filled  porous  media  is  presented  in  Chapter  II. 

The  implementation  of  these  theories  for  solution  on  a  digital  computer 
is  discussed  and,  finally,  the  modification  of  the  basic  normal-mode 
program  for  range-dependent  propagation  modeling  (i.e.,  adiabatic 
normal-mode  theory)  is  also  considered  in  Chapter  II.  The  normal-mode 
models  are  based  on  the  computer  model  of  Miller  and  Ingenito  (1975). 

The  experimental  procedures  are  discussed  in  Chapter  III.  The 
geometry,  type  of  experiments,  water  depth,  bottom  types,  etc., 
relating  to  each  of  the  experimental  at-sea  programs  is  first  discussed. 


A  description  of  the  data  collection  system  follows  and  includes  a 
brief  system  description,  bandwidths,  hyrophone  data,  etc.  The 
analysis  procedures  described  next  include  operations  performed  with 
the  Hewlett-Packard  Fourier  Analyzer  and  a  description  of  dispersion 
analysis  and  seismic-refraction  analysis.  Finally,  the  grab-sample  and 
core  analyses  performed  by  several  different  laboratories  are  discussed. 

The  analysis  of  experimental  data  begins  in  Chapter  IV.  The 
data  considered  in  this  chapter  is  from  a  test  site  off  the  coast  of 
Daytona  Beach,  Florida,  with  a  nominal  water  depth  of  19  m.  Also 
discussed  in  this  chapter  are  experimental  data  from  a  second  site  off 
Daytona  Beach  having  a  water  depth  of  30  m.  Experimental  data  from 
Corpus  Christi,  Texas,  for  a  water  depth  of  30  m  are  analyzed  in 
Chapter  V.  Chapter  VI  discusses  data  for  two  sites  with  nominal  water 
depths  of  18  and  30  m  off  the  coast  of  Panama  City,  Florida. 

The  analysis  in  Chapters  IV  through  VI  serves  to  establish  the 
validity  of  the  sediment  model  for  predicting  bottom  absorption.  In 
Chapter  VII,  transmission  loss  predictions  are  made  with  bottom 
absorption  obtained  from  the  sediment  models  for  a  region  off  the 
coast  of  Nova  Scotia  in  a  range-dependent  environment.  In  Chapter  VIII, 
the  study  is  summarized,  and  conclusions  are  made  concerning  the 
validity  and  failings  of  the  study. 


CHAPTER  II 


BACKGROUND 

2.1  Normal-Mode  Theory 

In  the  study  of  sound-transmission  problems,  the  solutions 
generally  fall  into  two  categories  based  primarily  upon  water  depth. 
Deep-water  problems  consider  those  situations  where  the  water  depth  is 
greater  than  approximately  100  fathoms.  Solutions  to  these  problems 
are  frequently  obtained  in  terms  of  rays  perpendicular  to  the  wavefronts 
which  emanate  from  the  source.  Shallow-water  problems  are  then 
concerned  with  water  depths  of  less  than  100  fathoms  and  frequently 
result  in  solutions  in  terms  of  the  normal  modes  set  up  in  the  water 
column  along  the  transmission  path.  As  might  be  expected,  the 
separation  between  the  two  types  of  solutions  is  not  that  clear;  for 
example,  very  low-frequency  waves  in  deep  water  may  yield  to  normal-mode 
solutions,  and  conversely,  high-frequency  waves  in  shallow  water  will 
obey  ray  theory.  Generally,  ray  theory  is  applied  where  the  wave  length 
is  a  small  fraction  of  the  water  depth.  For  the  water  depths  of 
interest  in  this  study  (less  than  33  fathoms),  both  types  of  solutions 
have  some  validity  since  the  frequencies  range  from  2  to  1000  Hz.  The 
emphasis  in  this  study  will  be  upon  normal-mode  solutions,  although  the 
ray  analogy  will  frequently  be  used  to  explain  normal-mode  theory. 

A  solution  in  terms  of  normal  modes  is  an  exact  solution  to  the 


wave  equation  for  certain  simple  boundary-value  problems.  If  the 
boundary  conditions  are  those  of  a  free  surface  and  a  rigid  bottom, 
we  have  perfect  reflection  at  the  surface  and  bottom,  and  an  ideal 
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waveguide  results.  The  solution  in  the  depth  coordinate  consists  of 
standing-wave  patterns  (e.g.,  as  in  an  organ  pipe),  or  the  normal  modes 
of  vibration  of  the  system.  For  shallow-water  propagation  in  the  ocean, 
the  problem  is  more  complex;  however,  the  same  concepts  apply.  The 
first  and  simpler  of  the  normal-mode  models  of  the  ocean  was  that  of 
Pekeris  (1948).  This  model  consists  of  an  isovelocity  water  layer  over 
another  fluid  layer  of  infinite  depth.  Densities  and  compressional 
sound  speeds  are  specified  in  each  layer,  and  the  water  depth  as  well 
as  the  other  parameters  must  be  invariant  with  range.  The  solution 
indicates  that  the  net  propagation  of  energy  is  along  the  channel  axis 
and  that  the  velocity  of  this  propagation  will  be  a  function  of 
frequency. 

To  elaborate  on  the  normal-mode  concept,  we  can  think  of 
shallow-water  propagation  in  terms  of  ray  theory.  The  singular 
characteristic  which  distinguishes  shallow-water  from  deep-water 
propagation  is  the  frequency  interaction  of  the  sound  field  with  the 
ocean  surface  and  bottom.  If  we  consider  propagation  at  some  distance 
from  the  source,  the  wave  fronts  will  strike  the  surface  and  bottom  as 
plane  waves.  For  each  encounter,  the  wave  will  undergo  a  phase  shift  at 
the  surface  and  bottom.  The  formation  of  standing  waves  requires  con¬ 
structive  interference  between  consecutive  downgoing  (or  upgoing)  rays 
striking  the  bottom  at  the  same  angle.  For  this  to  occur,  the  total 
phase  shift  due  to  surface,  bottom,  and  path  length  must  be  some 
multiple  of  360°.  Since  the  phase  shift  due  to  the  path  length  is  a 
function  of  frequency,  a  set  of  discrete  (harmonically  related) 
frequencies  will  satisfy  the  above  criteria  for  any  incidence  angle. 
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Other  frequencies  will  result  in  destructive  interference  for  this 
incidence  angle.  Thus,  the  standing-wave  patterns  set  up  in  the  sound 
channel  are  referred  to  as  the  normal  modes  of  the  channel.  The 
propagation  velocity  along  the  channel  axis,  termed  the  phase  velocity, 
is  dependent  upon  the  incidence  angle  through  a  simple  relation;  thus, 
a  particular  phase  velocity  is  associated  with  a  group  of  discrete 
frequencies.  The  multiple  reflections  and  consequent  interference  which 
develops  thus  causes  a  variation  in  propagation  velocity  with  frequency, 
termed  geometrical  dispersion.  With  a  similar  argument,  the  combined 
effects  of  all  the  preferred  angles  for  a  particular  frequency  may  be 
summed  to  create  the  mode-amplitude  distribution  in  the  water  column. 
That  is,  for  any  frequency,  the  pressure  amplitude  will  vary  with  depth 
in  a  predictable  manner,  this  variation  being  different  for  each  mode. 
The  only  loss  predicted  by  the  theory  is  inversely  proportional  to  the 
square  root  of  the  range,  termed  cylindrical  spreading. 


2.1.1  Formal  development  of  normal-mode  theory.  In  many  texts 
[Officer  (1958),  Tolstoy  and  Clay  (1966)],  the  formal  derivation  of  the 
normal -mode  solution  starts  with  the  general  wave  equation. 


V24> 


i 

C2  3t2 


for  the  velocity  potential  $ 
that  is. 


If  the  time  dependence  is  harmonic. 


<J>(t) 


-iwt 

e 


then  the  Helmholtz  equation  which  results  is  only  a  function  of  the 

space  variables,  giving: 

.  2 

V  <j>  +  ~  <j>  “  0 
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At  this  point,  many  derivations  assume  a  cylindrical  coordinate  system 
(r,  0,  2),  with  symmetry  in  the  0  coordinate,  and  assuming  range 
independent  parameters ,  the  separation  of  variable  approach  may  be  used 
to  obtain  equations  in  r  and  z  .  The  equation  in  r  is  Bessel's 
equation  which  leads  to  a  Bessel  function  solution  (for  outward 
propagation  of  energy)  for  the  r  coordinate.  The  z  solution  is 
that  of  an  ordinary  second-order  differential  equation  with  constant 
coefficients.  It  is  assumed  that  the  sound  velocity  c  is  a  function 
of  only  the  depth  coordinate.  The  boundary  conditions  for  the  Pekeris 
model  are  that  of  a  pressure  release  surface  at  the  air-water  interface, 
where  the  z  dependent  solution  must  go  to  zero ,  and  continuity  of 
pressure  and  the  vertical  velocity  component  across  the  water-bottom 
boundary.  Within  the  bottom  half-space,  the  solution  must  decay  to  zero 
with  Increasing  depth. 

The  source  conditions  are  applied  by  deriving  a  frequency- 
dependent  constant  which  matches  the  spherical  point  source  to  the 
complete  cylindrical  solution  of  the  wave  equation.  As  the  complete 
solution  involves  integration  of  the  product  solutions  for  r  and  z 
in  the  k  (wave  number)  plane  for  multi-valued  k  ,  complex  integration 
using  Cauchy's  residue  theorem  is  required.  The  residue  sum  for  this 
integration  in  the  complex  plane  then  yields  the  normal  modes. 

A  more  direct  derivation  of  the  normal -mode  equation  is  presented 
by  Ingenito  et  al.  (1978).  The  Pekeris  model  may  be  used,  but  to 
provide  a  better  introduction  for  the  computer  model,  the  three-layer 
normal-mode  model  in  Figure  1  will  be  assumed.  In  this  model,  all 
layers  are  assumed  t'  ’-e  fluid  layers;  hence,  compressional  velocities 
[c^(z),  C2(z),  c^l  are  assumed  throughout.  A  cartesian  coordinate 


A 
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system  will  be  used  for  most  of  the  derivation,  and  a  transformation 
to  cylindrical  coordinates  will  be  applied  near  the  end.  The  ocean 
surface  is  assumed  to  lie  in  the  x,y  plane,  while  the  z  coordinate 
increases  with  depth. 

Assume  a  unit  strength  harmonic  point  source  of  angular  frequency 
a)  at  a  depth  zq  as  shown.  If  the  time  dependence  is  dropped,  the 
velocity  potential  <}>(x,y,z)  must  satisfy  the  equation 

?  r  i2 

V‘ <fr  +  <j>  =  -6(x)6(y)6(z  -  zq)  .  (1) 

The  dimensional  requirements  of  Equation  (1)  may  be  satisfied  by  noting 

that  in  the  absence  of  boundaries ,  <p  will  have  the  solution 
1  j  If  R  /2 —  2  2 

<J>  =  e1  ,  where  R  =  /x  +  y  +  z  ,  and  both  sides  of  Equation  (1) 
will  have  the  same  dimensions.  We  may  then  find  $  at  any  point 
(x,y,z)  relative  to  the  velocity  potential  at  a  unit  distance. 

The  boundary  conditions  to  be  satisfied  are  those  stated  earlier. 
That  is,  at  the  surface, 

<j>(x,y,Q)  =  0  ,  (2) 

to  satisfy  the  pressure-release  condition.  At  the  boundary  between 
medium  1  and  2,  where  z  =  ,  we  have: 

PyJ)^  (x,y,H1)  =  (x.y,^)  (3) 


to  satisfy  the  condition  on  the  continuity  of  the  vertical  component 
of  velocity.  At  z  =  ,  in  a  similar  manner,  the  conditions  are: 
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and 


P24>(2)(x, 


Hx  +  H2)  = 


(3) 

P3<T  '(x,  y, 


Hx  +  H2) 


(5) 


36 


(2) 


3z 


z=H1+H2 


36 


(3) 


Z-Hl+H2 


(6) 


The  notations  ,  and  identify  the  velocity  potential 

in  the  three  regions  of  Figure  1.  The  two  final  conditions  are  that 
radiation  is  assumed  outward  from  the  source,  and  the  field  must  decay 
to  zero  in  the  half-space  of  medium  3. 

To  solve  Equation  (1),  the  double  Fourier  transform  will  be  used. 
That  is. 


and 
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u(nx,iyz) 
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<j>(x,y,z) 
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-i(n  x+n  v) 
e  y  dxdy  . 


(8) 


Substituting  Equation  (7)  into  Equation  (1)  yields  an  integro- 
differential  equation  which  may  be  solved  using  Equation  (8)  and  the 

shifting  property  of  the  Dirac  Delta  function.  Thus,  noting  that 

2  2  2  2  2  2 
'4f  =  n  +  n  *  and  k  (z)  =  u)  /c  (z)  ,  we  have: 

x  y 

2 

^4  +  [k2(z)  -  ip2]  U  =  -  5(z  -  z  )  ,  (9) 

dz  0 

a  non-homogeneous  differential  equation  for  U(nx>  f|y,  z)  .  To  solve 

Equation  (9),  it  is  assumed  that  Un(z)  is  the  eigenfunction  and  that 

k^  is  the  eigenvalue  of  the  following  equation: 


(10) 
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where  both  must  satisfy  the  boundary  conditions  (2)  through  (6). 

2  2  2 

Letting  \>n  *  k  (z)  -  ,  and  following  the  procedure  of  Clay  and 

Medwin  (1977),  we  find  that  for  the  nC^  and  m^  eigenfunctions  to 
satisfy  the  boundary  conditions  in  each  layer,  the  following  condition 
must  hold: 


I  pU  (z)U  (z)dz 
J  n  m 


6 

mn 


(ID 


o 

That  is,  the  eigenfunctions  Un(z)  form  an  orthonormal  set  satisfying 

Equation  (11),  where  p  will  assume  the  value  specified  in  each  layer. 

At  this  point,  the  meaning  of  and  its  effect  on  the 

eigenvalue  spectrum  (v  )  must  be  considered.  Physically,  k  =  —  , 

n  n 

where  v^  is  the  phase  velocity,  and  k^  is  the  wave  number  in  the 
direction  of  propagation.  The  spectrum  may  be  either  discrete  or 

continuous.  If  we  assume  that  c^  >  max  [c^,  c^]  ,  then  the  spectrum 
will  be  discrete  if 


c3  >  vn  >  max[c^(z) ,  c2(z)]  .  (12) 

The  discrete  spectrum  represents  energy  which  is  trapped  in  the  sound 
channel  due  to  perfect  reflection  at  the  surface  and  bottom.  If 
Equation  (12)  is  not  satisfied,  the  spectrum  will  be  continuous,  but 
these  modes  will  be  highly  damped  with  increasing  range  since  they  occur 
due  to  partial  reflection  at  the  bottom.  Thus,  the  continuous  modes  are 
important  at  only  very  short  ranges  and  will  hence  be  neglected  in  all 
further  discussions. 

Equation  (9)  may  now  be  solved  by  assuming  a  solution  composed  of 


an  expansion  of  the  eigenfunctions;  that  is, 


13 


U(VVz)  =  I  An(\,ny)Un(z)  *  (13) 

The  constant  A^r^.riy)  taay  he  evaluated  using  the  orthonormal 
relation  in  Equation  (11)  as: 
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The  solution  to  Equation  (9)  is  then: 


U  (z)U  (z  ) 

u(nx,ny,z)  =  p(*n)  l  n  n  ° 
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Thus,  transforming  back  to  the  space  domain, 


<f>(x,y,z) 
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In  order  to  solve  Equation  (16),  we  must  evaluate  the  integral 
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in  terms  of  polar  coordinates  (r,6)  .  This  is  done  in  Appendix  A 
and  yields  the  result  that 


I  =  iTr2H  (1)(k  r) 
o  n 


The  desired  relation  for  the  velocity-potential  solution  in  cylindrical 
coordinates  is  thus: 


f  P<*  >  1  U  (z  )U  (z)H  U)(k  r)  *  (17> 

4  o  L  non  o  n 
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In  Equation  (17),  U  (z  )  is  the  mode-amplitude  function  at  the 

no 

source.  U^(z)  is  t*ie  node-amplitude  function  at  the  receiver,  and 

H  ^ (k  r)  is  the  zero-order  Hankel  function  of  the  first  kind.  At 
o  n 

long  ranges  from  the  source,  the  Hankel  function  may  be  replaced  by  its 
asymptotic  approximation.  That  is, 
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The  velocity  potential  as  a  function  of  time  at  long  ranges  can  then  be 
expressed  as: 
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As  can  be  seen,  the  field  at  a  range  r  and  a  depth  z  diminishes  as 
1/i^r  and  propagates  with  a  phase  velocity  of  w/k^  .  A  broadband 
source  such  as  an  explosive  charge  will  contain  many  frequencies  with 
propagation  velocities  u^  (the  group  velocity)  related  to  the  phase 
velocity  by: 


n 


9m 

9k 


The 


pressure  field  is  related  to  the  velocity  potential  by  the  equation 


9<t> 

p  =  P 


so  that  the  variations  in  the  pressure  field  with  depth  due  to  the  source 
are  proportional  to  Un(zQ)  and  those  due  to  the  receiver  are 
proportional  to  Un(z)  • 
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The  computer  implementation  of  the  normal-mode  calculation  is 

based  on  solving  the  eigenvalue  problem  in  each  of  the  three  layers. 

This  requires  solving  the  differential  equation,  i.e..  Equation  (10), 

for  the  boundary  conditions  at  the  upper  and  lower  interfaces.  In 

reality,  the  differential  equations  are  replaced  by  difference 

equations  for  a  numerical  solution  and  the  sound-speed  profile  c(z) 

is  approximated  so  that  each  of  the  layers  is  divided  into  many  layers, 

each  having  a  constant  sound  speed.  The  solution  is  initiated  in  the 

third  layer  where  an  exact  solution  involving  only  one  unknown  constant 

(3) 

may  be  written  (from  the  requirement  that  U  (z)  -*•  0  as  z  -*•<*>] . 

n 

(3) 

A  trial  value  of  k  is  assumed,  and  U  and  its  derivatives  are 

n  n 

calculated.  The  boundary  conditions  are  then  satisfied, and  the  solution 
proceeds  upward  within  each  layer  until  the  surface  is  reached  where 
Un(1) (0)  must  be  zero  or  less  than  some  pre-selected  small  number.  The 
procedure  is  repeated  until  this  last  condition  is  satisfied. 

The  transmission  loss  may  be  calculated  from  the  solution  obtained 
for  and  kn  .  That  is,  the  modal  pressure  sum  may  be  obtained  from 

Equation  (18),  and  then  the  energy  in  each  mode  calculated.  The 
incoherent  transmission  loss  expressed  in  decibels  for  unit  source 
strength  may  then  be  expressed  as : 
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10  log 


10 


2^p2(zQ) 


U  (z  )U  (z) 
non 

A  r 


(19) 


where  U  (z) 
n 


and  k^  are  evaluated  numerically  at  the  receive  point 
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of  interest. 
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Certain  modeling  situations  involving  near-surface,  consolidated 
(and  maybe  some  semi-consolidated)  bottom  structures  with  high  shear 
velocities  require  a  solid-bottom  model  as  opposed  to  the  fluid  model 
just  described.  Ingenito  et  al.  (1978)  derives  this  type  of  model, 
replacing  the  fluid  half-space  (layer  3)  with  a  solid  half-space.  The 
primary  difference  in  the  solution  is  that  two  eigenfunctions  must  now 
be  determined  in  the  half -space,  one  for  the  compressional  wave  and  one 
for  the  shear  wave.  The  boundary  conditions  at  z  *  also 

involve  the  shear  eigenfunction  and  must  force  this  solution  to  become 
zero  within  the  upper  two  fluid  layers.  A  complete  description  of  the 
implementation  of  both  fluid  and  solid  normal-mode  computer  models  is 
presented  by  Miller  and  Ingenito  (1975). 

2.1.2  Loss  mechanisms.  In  the  preceding  discussion,  the 
assumption  was  made  that  for  large  r  ,  all  of  the  energy  was  trapped 
within  the  sound  channel.  This  resulted  from  the  assumption  of  perfect 
reflection  at  the  water-sediment  interface.  The  model  thus  predicts 
transmission  loss  due  only  to  geometrical  spreading  (I.e.,  l//r  or 
cylindrical  spreading).  As  frequently  happens,  however,  this  model  of 
the  ocean  and  bottom  is  too  simple.  A  physically  realistic  ocean 
bottom  will  absorb  energy  from  all  the  Incident  rays;  shear  waves  will 
be  generated  due  to  the  incident  compressional  waves;  and  finally, 
bottom  roughness  will  cause  energy  to  be  scattered  at  unfavorable 
angles  resulting  in  attenuation.  The  net  effect  is  that,  due  to  these 
mechanisms,  propagation  losses  in  the  real  ocean  environment  are  much 
greater  and  have  more  variability  than  the  basic  theory  predicts. 


As  noted,  three  loss  mechanisms  associated  with  the  bottom 
have  been  identified;  that  is,  absorption  of  compressional  waves, 
scattering  losses,  and  conversion  of  compressional  waves  to  shear 
waves.  Losses  due  to  bottom  absorption  were  noted  in  Chapter  I;  and, 
since  they  are  of  primary  interest  in  this  study,  they  will  be  discussed 
at  length  near  the  end  of  this  chapter.  The  other  loss  mechanisms  are 
described  briefly  in  the  following  paragraphs. 

Loss  due  to  scattering  can  be  viewed  as  resulting  from  the 
transfer  of  energy  between  propagating  normal  modes  due  to  the  roughness 
of  the  ocean  floor.  McDaniel  (1977)  has  analyzed  this  loss  mechanism 
and  shows  that,  for  the  higher  frequencies  (greater  than  200  Hz),  the 
loss  due  to  mode  conversion  can  be  of  the  same  magnitude  as  volume 
absorption.  The  loss  results  from  the  transfer  of  energy  to  higher- 
order  modes,  which  in  turn  have  larger  attenuation  coefficients 
resulting  in  greater  loss  than  if  conversion  were  not  present.  A 
sloping  bottom  may  also  cause  this  same  type  of  loss  (i.e.,  mode 
conversion);  however.  Pierce  (1965)  has  shown  that  if  the  water  depth 
varies  slowly  with  range,  mode  conversion  may  be  neglected.  The 
basic  normal-mode  theory  may  be  modified  to  handle  this  situation. 

The  modified  theory  is  termed  adiabatic  normal-mode  theory  and  is 
discussed  in  the  latter  part  of  this  chapter. 

Loss  due  to  conversion  to  shear  waves  will  be  present  for  those 
sediments  capable  of  supporting  shear  waves.  Hamilton  (1969)  has  found 
that  most  porous,  marine  sediments  will  support  shear  waves.  The  loss 
due  to  conversion  to  shear  waves  is  treated  by  Officer  (1958),  Eby  et 
al.  (1960),  and  Williams  and  Eby  (1962)  by  deriving  a  plane-wave 
reflection  coefficient  for  the  case  of  a  rigid  bottom,  and  then  using 
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this  coefficient  to  include  a  loss  per  bounce  in  the  Pekeris  model.  For 

shear  velocities  on  the  order  of  400  m/sec  or  greater,  this  loss 

mechanism  may  be  the  largest  of  the  three. 

In  the  model  for  normal-mode  propagation,  all  layers  are 

considered  lossless.  In  order  to  account  for  the  attenuation  associated 

with  real  ocean  sediments,  the  bottom  layers  must  be  modified  to 

represent  a  lossy  medium.  This  is  done,  following  Kornhauser  and 

Raney  (1955),  and  Tolstoy  (1958),  by  multiplying  the  basic  solution 

-6nr 

[Equation  (18)]  by  an  exponential  attenuation  factor,  e  ,  where  6^ 
is  the  mode-attenuation  coefficient  .  The  modal  sum  for  the  velocity 
potential  then  becomes: 
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It  may  be  shown  [Kornhauser  and  Raney  (1955),  and  Ingenito  (1973)]  that 
6^  can  be  equated.  If  we  assume  all  the  attenuation  is  due  to  the 
absorbing  bottom,  to  the  product  of  a  bottom-absorption  coefficient  a 

A 

and  a  modal  energy  ratio,  y  ,  equal  to  the  fraction  of  the  total  energy 

n 

(2) 

trapped  in  the  bottom.  That  is,  for  the  second  layer,  6  *  Y  a., 
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With  this  background,  the  scope  of  the  study  can  be  narrowed  to 
specify  one  aspect  of  the  problem.  That  is,  if  the  mode-attenuation 
coefficient  for  any  mode  (e.g.,  the  first  mode)  is  known,  and  if  a 
computer  implementation  of  the  normal-mode  equations  is  available  so 
that  may  be  calculated,  then  the  bottom-absorption  coefficient  may 

be  calculated  as  well. 

In  general,  the  experimental  data  obtained  at  sea  will  yield 
transmission  loss  as  a  function  of  range.  To  generate  this  information, 
received  energies  are  calculated  at  a  particular  frequency  and 
subtracted  from  known  source  levels  at  each  range  point  of  interest, 
the  difference  being  the  transmission  loss.  The  mode-attenuation 
coefficient  may  be  calculated  for  any  mode  and  frequency  by 
(1)  first  isolating  the  received  energy  in  that  mode,  (2)  calculating 
the  transmission  loss  for  the  mode,  (3)  applying  a  correction  for 
cylindrical  spreading  losses,  and  (A)  computing  the  slope  of  the 
corrected  transmission  loss  versus  range  curve.  The  computed  slope, 
with  units  of  dB/km,  is  the  mode-attenuation  coefficient  for  the 
isolated  mode  and  the  frequency  of  interest.  A  summary  of  the 
research  documenting  the  use  of  these  techniques  is  presented  in 
the  following  paragraphs. 

2.1.3  Previous  shallow-water  studies.  Although  the  study  of 
mode  attenuation  or  bottom  loss  was  not  the  primary  concern,  the  work 
of  Worzel  and  Ewing  (19A8)  was  probably  the  first  experimental  work  at 
sea  with  explosive  charges  to  be  tested  against  normal-mode  theory. 
During  the  same  period,  Pekeris  (1948)  had  developed  normal-mode  theory 
for  sound-propagation  problems  in  shallow  water  and  was  able  to  explain 
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many  of  Che  observed  results  of  Worzel  and  Ewing  with  the  theory. 

Eby  et  al.  (1960)  performed  model  experiments  with  a  water  layer  over 
a  Hycar  rubber  bottom.  Over  the  range  of  test  frequencies,  55  to 
600  kHz,  experimental  mode-attenuation  coefficients  compared  well  with 
theory  at  the  low  end  of  the  range,  but  toward  the  upper  end,  losses 
attributed  to  conversion  to  shear  waves  caused  higher  than  expected 
attenuation.  Tolstoy  (1958)  performed  normal-mode  experiments  in 
22.6  m  of  water  off  the  Atlantic  coast  and  found  that  mode-attenuation 
coefficients  for  the  first  mode  ranged  from  3.57  to  5.13  dB/km  over  a 
frequency  range  of  88  to  148  Hz.  Kibblewhite  and  Denham  (1966) 
performed  experiments  in  water  depths  ranging  from  91  to  274  m  off 
the  New  Zealand  coast.  Their  results  led  them  to  conclude  that  the 
attenuation  was  dominated  by  absorption  and  conversion  to  shear  at  low 
frequencies  (less  than  50  Hz),  and  that  scattering  dominated  around 
100  Hz  (where  an  attenuation  minimum  occurred)  and  also  at  higher 
frequencies.  More  recently,  Ferris  (1972)  and  Ingenito  (1973) 
collaborated  on  measurements  performed  in  the  Gulf  of  Mexico  at 
frequencies  of  400  and  750  Hz.  They  obtained  good  agreement  between 
theoretical  and  experimental  mode-amplitude  distributions  for  the  first 
two  modes  along  a  constant-depth  (31  m)  propagation  path  over  a  hard- 
sand  bottom.  Finally,  Ingenito  and  Wolf  (1976)  made  measurements  off 
the  Yucatan  Peninsula  along  a  constant-depth  (30  m)  track  overlying  a 
near-surface  limestone  layer.  Conversion  due  to  shear  was  found  to 
be  the  dominant  loss  mechanism  for  this  high-velocity  layer  (1900  m/sec 
compressional  and  1000  m/sec  shear). 
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2.1.4  Range-dependent  modeling.  A  basic  assumption  involved 
in  the  derivation  of  the  normal-mode  equation  [Equation  (18)  ]  was  that 
all  parameters  were  independent  of  range,  permitting  the  separation  of 
range  and  depth  solutions.  While  this  assumption  is  not  unduly 
restrictive,  there  being  many  situations  where  sound-velocity  profiles, 
water  depths,  etc.,  do  not  vary  substantially  with  range,  it  is  often 
desirable  to  study  propagation  under  range-dependent  conditions.  To 
apply  normal-mode  theory  to  these  situations,  the  adiabatic  approxima¬ 
tion  is  used;  that  is,  it  is  assumed  that  each  normal  mode  propagates 
independently  with  no  exchange  of  energy  between  modes  (i.e. ,  no  mode 
coupling).  As  noted  earlier,  McDaniel  (1977)  found  that  losses  due  to 
mode  coupling  increased  for  the  higher  frequencies;  hence,  it  is 
expected  that  use  of  the  adiabatic  approximation  will  be  limited  to 
the  lower  frequencies.  The  range-dependent  calculation  of  transmission 
loss  is  then  made  using  the  range-independent  model  with  the  assumption 
that  the  modes  adapt  to  the  environment  at  the  point  of  interest,  i.e., 
locally,  propagation  is  independent  of  range. 

These  assumptions  have  been  tested  in  experimental  work  by 

Eby  et  al.  (1960),  and  Ingenito  et  al.  (1978),  and  treated  theoretically 

by  Pierce  (1965).  The  range-dependent  model  may  be  implemented 

(Ingenito  et  al.)  by  modifying  Equation  (20)  so  that  the  modal  wave 

number  k  and  attenuation  coefficient  5  are  replaced  bv  comparable 
n  n 

quantities  k  and  6  which  are  averaged  over  the  propagation  path. 
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and 
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The  expression  for  the  incoherent  transmission  loss  then  becomes 
[following  Equation  (19)]: 
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The  procedure  is  then  to  specify  parameters  for  the  water  and  seabed 
at  a  number  of  range  points  along  the  propagation  track.  The  range- 
independent  model  discussed  previously  is  run  at  these  range  points  and 
U  (z  ),  k  ,  and  5  calculated.  The  receive  function  U  (z)  is 
calculated  only  at  the  zero  range  point.  If  ranges  are  required  other 
than  those  at  which  the  water  and  bottom  parameters  have  been  specified. 
Interpolation  of  ^n(z0)  between  specified  range  points  is  required. 


2.2  Geoacoustic  Models  of  the  Ocean  Bottom 

In  specifying  models  of  the  ocean  bottom,  the  principal  acoustic 
parameters  of  interest  are  the  density,  compressional  and  shear 
velocities,  and  the  bottom-absorption  coefficient.  Related  to  these 
acoustic  properties  of  the  model,  by  largely  empirical  formulations, 
are  engineering  properties  such  as  porosity,  grain-size  and  distribu¬ 
tion,  and  the  elastic  constants  (such  as  the  bulk  modulus  and  rigidity 
modulus).  In  addition,  a  complete  geoacoustic  model  details  the 
thickness  and  properties  of  the  sediment  and  rock  layers  down  to  a 
depth  which  is  dependent  on  the  lowest  frequency  of  interest. 
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Several  different  techniques  are  required  to  formulate  such  a 
model.  For  the  deeper  layers,  the  required  acoustic  parameters  must  be 
determined  by  indirect  methods  such  as  seismic-refraction  and 
dispersion  analysis,  reflection  profiling,  or  literature  surveys.  The 
surficial  sediments  may  be  sampled  through  the  use. of  corers  and  grab 
samples.  Analysis  of  the  cores  or  grab  samples  will  provide  many  of  the 
engineering  properties  such  as  grain-size  data  (and  porosity  for  cores) . 
With  this  data,  absorption  coefficients  may  be  predicted  using  a 
sediment  model  for  the  upper  few  meters  of  the  seabed. 

2.2.1  Sediment  model  of  Biot.  As  discussed  in  Chapter  I,  the 
work  of  Biot  (and  more  recently,  Stoll)  and  Hamilton  has  dominated  the 
literature  with  regard  to  sound  absorption  in  marine  sediments.  In  the 
following  paragraphs,  the  model  of  Biot  will  be  derived,  using  Stoll's 
development,  and  the  work  of  Hamilton  discussed.  In  Chapter  III,  the 
experimental  techniques  (such  as  seismic-refraction  and  dispersion 
analysis)  for  determining  sub-bottom  parameters  will  be  discussed. 

The  theory  developed  by  Biot  is  a  comprehensive  description  of 
the  response  of  linear,  porous  materials  containing  compressible  fluid. 
Stoll  clarifies  and  presents  an  abbreviated  version  of  Biot's  develop¬ 
ment  forming  the  basis  for  the  following  discussion. 

The  development  of  Stoll  begins  with  a  pair  of  coupled 
differential  equations  derived  by  Biot  for  the  motion  of  the  skeletal 
frame  and  the  motion  of  the  viscous  pore  fluid  relative  to  the  frame. 
That  is,  for  the  frame, 

V2(He  -  CO  =  —  j  (pe  -  pfO 
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and  for  the  fluid  relative  to  the  frame. 
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Solutions  for  e  ,  the  dilitational  strain  of  the  frame,  and  £  ,  the 
increment  of  fluid  flow,  are  assumed  to  be  of  the  form 
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(25) 


and 
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If  Equations  (25)  and  (26)  are  substituted  in  Equations  (23)  and  (24) , 
two  equations  in  and  result.  If  solutions  exist  for  A^  and 

A^  ,  it  must  hold  that  the  determinant  of  the  coefficients  of  A^  and 
A^  equals  zero.  Stoll's  determinant  relation  then  results  as: 
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A  solution  for  the  complex  variable  £  will  yield  the  absorption  (the 
imaginary  part)  and  phase  velocity  (the  real  part)  for  two  types  of 
compressions!  waves.  A  relation  similar  to  Equation  (27)  will  yield  a 
comparable  solution  for  the  shear  wave.  The  first  compressional  wave 
is  the  usual  one  associated  with  fluids;  the  second  wave  has  lower 
velocity  and  is  greatly  attenuated.  Recently,  Plona  (1980)  observed 
the  second  wave  at  ultrasonic  frequencies  in  a  fluid-saturated  porous 
medium  consisting  of  water  and  sintered  glass  spheres.  Only  the  first 
compressional  wave  is  considered  here. 


25 


If  the  determinant  in  Equation  (27)  is  reduced,  the  following 
fourth-order  complex  equation  is  obtained: 


a/4  +  a2*2  +  a3  = 

The  constants  a^  ,  ,  and  a^ 

follows: 


0  .  (28) 
are  complex  and  are  defined  as 


ax  =  C2  -  H  M  ,  (29) 

a2  =  (Hm  +  pM  -  2Cpf)u>2  -  i  ^  H  (30) 

and 

a3  «  (Pf2  -  pm)(i/*  +  i  p  .  (31) 

In  Equations  (29),  (30),  and  (31),  p  is  the  saturated  density;  p^  is 
the  density  of  the  pore  fluid;  K  is  the  sediment  permeability;  0  is 
the  fluid  viscosity;  and  <±)  is  the  radian  frequency.  The  complex 
constants  (C  ,  H  ,  and  M)  are  functions  of  fluid,  frame,  and  frame 
material  bulk  moduli  (k^  ,  k^  ,  and  k^  ,  respectively)  frame  shear 
bulk  modulus  U  ,  porosity  8  ,  and  a  parameter  m  representing  the 
tortuosity  of  the  pore  spaces.  The  constants  are  defined  by  Stoll  as 
follows: 
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In  Equations  (30)  and  (31) ,  the  term  represents  the  frequency- 

dependent  viscous  effects  through  the  correction  factor  F  derived  by 


Biot.  The  correction  factor  F  is  defined  as  a  function  of  <  ,  where 
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The  functions  ber(<)  and  bei(K)  are  real  and  imaginary  parts  of  the 
Kelvin  function,  and  a  is  a  pore-size  parameter.  The  functions 
ber'(K)  and  bei' (K)  are  the  derivatives  of  the  Kelvin  functions. 

To  account  for  the  unknown  nature  of  the  pore  size  and  shape, 
the  parameters  m  and  a  are  introduced.  The  parameter  m  is  defined 
as : 

m  -  apf/8  , 

where  a  equals  one  for  a  system  of  uniform  pores  with  axes  parallel  to 
the  pressure  gradient  and  three  for  a  random  orientation  of  uniform 
pores.  In  practice,  a  must  be  treated  as  an  experimental  parameter. 
The  parameter  a  is  related  to  pore  size  and  permeability  and  will  be 
discussed  later.  The  essential  difference  between  a  and  m  is  in 
their  usage;  m  has  an  inertial  effect  on  flu:H  flow,  reducing  this 
flow  for  the  more  random  pore  orientation;  a  affects  the  viscous 
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resistance  to  fluid  flow  and  results  in  an  increased  attenuation  for 
high  frequencies. 

The  inelastic  nature  of  the  frame  is  included  in  the  model  by 
making  the  frame  bulk  modulus  k^  and  the  frame  shear  modulus  y 
complex: 


The  Young's  modulus  for  longitudinal  vibration  (E)  is  calculated  from 
experiments  on  the  longitudinal  vibration  of  a  thin  rod  of  the  frame 
material  (also,  see  Hardin  and  Richart,  and  Hardin),  where  E  governs 


the  longitudinal  velocity  of  propagation  v£  through  the  relation: 


v 


E 


(43) 


The  complex  part  E*  is  obtained,  in  a  manner  similar  to  y*  ,  from 
the  log  decrement  for  longitudinal  vibrations;  that  is, 


A 


E 


7TE* 

E 


(44) 


Thus,  in  summary,  k^  and  y  are  obtained  from  experimental  values  of 
As  ’  AE  *  vs  ’  and  VE  * 

It  is  necessary  to  use  several  data  sources  in  order  to  obtain 
the  remaining  input  parameters.  The  bulk  moduli  of  the  fluid  and  frame 
material  can  be  approximated  from  handbook  data  once  the  frame  material 
is  determined  by  core  analysis.  The  porosity  may  be  estimated  from  data 
in  Hamilton  (1972),  where  porosity  is  plotted  as  a  function  of  mean 
grain-size  if  data  from  core  analysis  is  unavailable.  The  fluid  and 
frame  material  density  may  be  obtained  from  a  handbook,  and  the  saturated 
density  for  the  medium  may  then  be  calculated  from  the  equation 


P  =  (1  -  9)Pr  +  3pf 


(45) 


The  absolute  viscosity  of  the  pore  fluid  n  is  a  standard  value 
2 

(0.01  dyne-sec/cm  at  20°C)  obtainable  from  a  handbook. 

In  order  to  run  the  sediment  model,  k^_  ,  k^  ,  ,  P  ,  K,  *  , 

a  ,  v  ,  v  ,  y  ,  0  ,  A  ,  and  A_  must  be  supplied  as  input  data, 
c<  o  ti  S 

along  with  the  desired  frequency.  The  permeability  K  iu  calculated 
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from  the  sediment  mean  grain-size  (d  in  mm)  and  phi  deviation 
measures  (o^  in  phi  units)  through  the  equation 

K  =»  7.6d2e“1,31°0  x  10~6  (46) 


from  Krumbein  and  Monk  (1942).  The  pore-size  parameter  a  is  then 
calculated  from  the  relation 


a  -  3.87  v'K 


(47) 


used  by  Stoll.  The  bulk  modulus  of  the  frame  material  k  mav  be 

r 

calculated  using  the  Voigt-Reuss-Hill  (VRH)  averaging  method  described 
by  Hamilton  (1969).  That  is,  if  the  frame  is  composed  of  several 
mineral  types,  an  average  value  based  on  the  fractional  (by  volume) 
concentration  of  each  may  be  calculated  from  the  equation  (e.g.,  for  a 
two- component  frame) : 


kLk2(l  +  +  V22)  +  V1V2(k12  +  k22) 

2k2Vl  +  2klV2 


(48) 


where  k^  and  k2  are  the  bulk  moduli  of  the  two  frame  materials,  and 
and  V2  are  the  fractional  concentrations  of  each. 

Finally,  once  the  above  parameters  have  been  specified  as  input 
to  the  main  program,  computer  subroutines  are  used  to  calculate  the 
Kelvin  functions,  their  derivatives  and  the  roots  of  Equation  (28). 


2.2.2  Hamilton's  approach  to  sediment  absorption.  The  work  of 
Hamilton,  while  mostly  empirical,  is  based  to  some  extent  on  linear 
viscoelastic  theory  as  presented  by  Ferry  (1961).  The  linear  visco¬ 
elastic  theory  assumes  that  a  complex  ratio  between  stress  and  strain 
exists,  and  the  damping  or  absorption  is  then  represented  by  the  phase 
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angle  of  the  ratio.  The  absorption  is  predicted  to  vary  as  the  first 
power  of  frequency. 

The  empirical  work  of  Hamilton  derives  from  many  measurements 
on  cores,  artificial  sediments,  and  short-range  measurements.  The 
results  cover  a  frequency  range  of  5  Hz  to  1  MHz,  although  the  greatest 
amount  of  data  is  above  1  kHz.  The  work  is  described  in  Hamilton 
(1972,  1974,  1975).  The  technique  for  predicting  the  absorption  of 
compressional  waves  is  best  described  in  Hamilton  (1972) ,  where 
regression  equations  for  absorption  as  a  function  of  mean  grain-size 
and  porosity  are  presented.  The  absorption  data  is  presented  as 
(*A  =  K^f^  ,  with  Ka  having  units  of  dB/ (m  •  kHz)  and  f  expressed  in 
kHz  so  that  is  in  dB/m.  In  general,  the  absorption  coefficient 

aA  will  be  expressed  in  units  of  dB/km.  As  an  example  of  the  data 
given  by  Hamilton,  Table  1  is  extracted  from  the  last  reference.  Given 
any  mean  grain-size  and  frequency,  the  absorption  coefficient  may  be 
predicted  according  to  Hamilton  through  the  data  in  Table  1. 

Data  exists  which  appears  to  support  both  theories,  the  source  of 
this  data  being  measurements  in  cores,  artificial  sediments  and  in-situ 
measurements  over  very  short  ranges  (1  to  2  m) .  The  absorption 
measurements  on  cores  and  artificial  sediments  can  be  considered  as 
measures  of  the  intrinsic  absorption  in  the  sediment  and  rock  material, 
while  the  in-situ  measurements  are  considered  as  the  "effective 
attenuation"  measured  from  all  sources. 

Hamilton  (1972  and  1975)  has  accumulated  a  very  large  body  of 
data  over  a  frequency  range  of  5  Hz  to  1  MHz  which  appears  to  show 
nearly  a  linear  dependence  on  frequency.  This  data  includes  both 
intrinsic  measurements  [McCann  and  McCann  (1969),  for  example]  and 
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TABLE  1 

ABSORPTION  REGRESSION  EQUATIONS 
[Hamilton  (1972)] 

Mz  Range 
(phi*) 

0  -  2.6 
2.6  -  4.5 
4.5  -  6.0 
6.0  -  9.5 


Regression  Equation 
for  R. 

_ A _ 

Ka  -  0.4556  +  0.0245  M 
A  z 

K.  -  0.1978  +  0.1245  M 
A  z 

Ka  -  8.0399  -  2.5228  M 
A  z 

Ka  =  0.9431  -  0.2041  M 
A  z 


* 

The  phi  unit  is  defined  as  M  ■  -  log  d  , 

z  z 

where  d  is  the  mean  grain  diameter  in  mm. 


+  0.20098  M  2 
z 

-  0.0117  M  2 
z 
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effective  attenuation  measurements  [Neprochnov  (1971)].  Acoustic 
experiments  conducted  at  sea  by  MacKenzie  (1960)  and  Cole  (1965) 
also  appear  to  support  Hamilton's  theory. 

Literature  supporting  the  work  of  Biot  and  Stoll  includes 
laboratory  experiments  by  Busby  and  Richardson  (1957),  Shumway  (I960), 
Nolle  et  al.  (1963),  Hampton  (1967),  and  Johnston  et  al.  (1979).  The 
work  of  Hampton  also  reports  measured  velocity  dispersion  of  approx¬ 
imately  two  percent  over  the  frequency  range  of  5  to  200  kHz.  Acoustic 
measurements  reported  by  Ingenito  (1973)  indicate  that  attenuation 
varies  as  with  frequency. 


CHAPTER  III 


EXPERIMENTAL  PROCEDURES 


3.1  Background 

In  order  to  study  the  effects  of  sediment  absorption  on 
transmission  loss,  two  types  of  information  must  be  obtained 
experimentally.  First  of  all,  to  measure  the  effectiveness  of 
techniques  developed  for  modeling  the  sediment  and  predicting  its 
absorption,  experimental  sediment-absorption  data  must  be  available 
for  the  frequencies  of  interest.  As  discussed  briefly  in  Chapter  I, 
this  experimental  data  is  obtained  by  measuring  the  transmission  loss 
in  regions  having  a  wide  variety  of  bottom  types.  From  the  transmission- 
loss  measurements,  the  attenuation  attributable  to  sediment  absorption 
may  be  assessed.  Thus,  the  test  sites  and  the  methods  used  to  obtain, 
process,  and  analyze  the  experimental  data  are  of  great  importance. 
Secondly,  the  surficial  sediments  and  sub-bottom  structure  must  be  well 
known  so  that  sediment  and  sub-bottom  modeling  efforts  will  represent 
the  bottom  to  the  highest  degree  possible.  Thus,  good  grab-sample  and 
coring  data  are  required  to  characterize  the  surficial  sediments. 
Experiments,  such  as  the  seismic-refraction  experiment,  have  been 
designed  so  that  the  data  obtained  will  yield  information  on  the 
sub-bottom  structure.  Analysis  of  the  dispersive  propagation  effects 
provides  another  technique  for  characterizing  the  bottom  parameters. 
Finally,  reflection  profiling  may  provide  both  surficial  sediment  data 


and  evidence  as  to  the  shallow  sub-bottom  structure. 
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With  this  in  mind,  test  sites  were  selected  to  provide  a 
variety  of  sediment  types;  in  most  cases,  regions  were  selected  with 
constant  depth  and  sloping  bottoms  having  nominal  water  depths  of  15  to 
46  m.  The  Scotian  Shelf  test  area  was  the  one  exception  to  these 
criteria,  as  a  sloping  bottom  in  deeper  water  was  desired.  The  test 
sites  selected  are  listed  in  Table  2  along  with  their  general  bottom 
type.  As  can  be  seen,  a  variety  of  bottom  types  are  present. 

In  the  following  sections,  the  experimental  geometry,  sound 
sources,  instrumentation,  analysis  system  and  procedures,  and  the  source 
and  types  of  sediment  data  used  are  discussed.  This  information,  along 
with  the  theory  discussed  in  Chapter  II,  provides  the  necessary 
background  for  understanding  the  analysis  and  comparisons  of  data  in 
the  later  chapters. 

3.2  At-Sea  Experiments 

From  1973  through  1978,  experimental  programs  were  conducted  to 
measure  the  characteristics  of  acoustic  waves  propagating  in  shallow 
water.  In  addition,  environmental  measurements  were  made  to  determine 
seabed  properties.  Tests  were  performed  at  two  different  sites  off 
Daytona  Beach  during  the  1976  tests  and  at  one  site  during  the  1973 
tests.  Data  is  also  reported  for  two  sites  off  Panama  City,  one  on 
the  Scotian  Shelf,  and  one  off  Corpus  Christi. 

The  experiments  were  generally  of  three  different  types,  all  of 
which  involved  propagation  over  a  track  of  constant  depth  and  over  a 
track  from  a  source  in  deep  water  to  a  shallow-water  receiver. 
Seismic-refraction  experiments  were  designed  to  yield  information  about 


the  sub-bottom  structure.  Short-range  experiments  were  designed  to 
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TABLE  2 

TEST  SITES  AND  BOTTOM  TYPES 


Test  Site 

Panama  City,  Florida 
Daytona  Beach,  Florida 

Corpus  Christi,  Texas 
Scotian  Shelf 


Date  of  Tests 

June  1976 

November  1973 
November  1976 

April  1977 

September  1978 


Bottom  Type 
Sand 
Sand 

Mud 

Gravel 
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yield  mode-amplitude  and  mode-attenuation  data.  The  long-range 
experiments  permitted  calculation  of  mode-attenuation  data  at  longer 
ranges,  which  has  certain  advantages  if  only  the  first  mode  is  desired. 
Dispersion  analysis  was  performed  upon  selected  signals  from  all  three 
types  of  experiments.  The  experiments  are  described  as  follows: 

1.  Short-range  experiments  employed  small  7.1  g  (1/4  oz) 
charges  detonated  at  depths  selected  to  suppress  certain 
modes,  while  enhancing  others.  The  receiver  for  these 
tests  was  a  bottom-moored  vertical  array  of  hydrophones 
spanning  the  water  column  in  31  m  of  water.  The  maximum 
range  for  these  tests  was  limited  by  line-of-sight  for 
radio-frequency  propagation,  typically  12  m  for  the 
antenna  placement  used.  Shots  were  deployed  at  range 
intervals  of  approximately  1  km.  For  these  tests,  the 
shot  boat  also  served  as  the  recording  platform  for  data 
telemetered  from  the  receiver. 

2.  Experiments  to  study  propagation  at  longer  ranges  employed 
an  aircraft  (or  another  ship)  to  drop  larger  charges.  The 
aircraft  operations  typically  used  816  g  (1.8  lb  SUS) 
charges  detonated  at  depths  of  18.3  and  244  m.  The 
maximum  range  was  approximately  100  km,  and  the  shot 
spacing  in  range  was  nominally  18  km.  The  receiver  for 
these  tests  was  usually  a  single  hydrophone,  one  meter 
above  the  bottom,  the  output  of  which  was  telemetered  to 
the  support  ship  for  recording. 
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3.  The  seismic-refraction  experiments,  designed  to  yield 
sub-bottom  information,  used  454  g  (1  lb)  charges 
detonated  on  the  bottom.  The  maximum  range  and  spacing 
were  very  similar  to  the  short-range  tests.  The 
receiver  was  the  single  hydrophone  and  a  three-axis 
geophone;  again,  the  data  were  recorded  aboard  ship  as 
the  charges  were  dropped  at  increasing  ranges.  At  some 
sites,  the  experimental  procedure  varied  somewhat  from 
those  described  here — these  will  be  noted  as  they  occur. 

In  addition  to  the  acoustic  information  obtained  during  the 
experiments,  other  data  such  as  sound-velocity  profiles,  cores,  grab 
samples,  sea  conditions,  and  sub-bottom  profiles  were  obtained  whenever 
possible.  On  those  occasions  when  this  data  could  not  be  obtained 
during  the  experiment,  the  best  available  alternate  data  sources  were 
used. 

3.3  Data  Collection  System 

The  basic  requirement  of  the  data-collection  system  was  to 
receive  the  acoustic  signals  which  propagated  through  the  water  from  the 
source,  process  the  signals  so  they  were  suitable  for  recording  on 
analog  tape,  and  finally,  record  the  signals.  Two  basic  systems  have 
been  used  for  accomplishing  this:  (1)  a  system  using  a  telemetry  link 
to  transfer  data  from  the  sensors  to  the  recorder  aboard  a  ship,  and 
(2)  a  system  using  a  cable,  hardwiring  the  sensors  to  the  recorder 
mounted  on  an  oceanographic  platform.  The  first  system  was  used  for  all 
experiments  except  the  Panama  City  tests  where  the  second  system  was 
used.  Both  systems  will  be  discussed,  and  the  sensors  and  sources  used 


will  be  described. 
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3.3.1  Shipboard  and  in-water  instrumentation.  The  buoy  system 
is  pictured  in  Figure  2  and  consists  of  a  submerged  unit,  sensors,  an 
underwater  RF  cable,  and  a  small  antenna  buoy.  The  sensor  outputs  were 
processed  by  amplifying  the  signals,  with  preamplifiers  and  variable- 
gain  amplifiers,  and  then  filtering  the  signals  to  a  frequency  band 
dictated  by  the  telemetry  system.  The  processed  signals  were  then  used 
to  frequency  modulate,  in  an  FM/FM  process,  a  data  transmitter.  The 
transmitter  was  coupled  to  the  antenna,  mounted  on  a  surface  buoy,  by 
an  underwater  RF  cable.  The  submerged  unit  also  contained  batteries 
and  several  other  sensors  such  as  a  leak  detector,  battery  monitor,  and 
a  two-axis  pitch  sensor.  Data  were  received  at  the  ship  and  recorded 
on  analog  tape.  Commands  were  sent  from  the  ship  to  the  buoy  to  change 
the  amplifier  gains  and  other  buoy  functions.  System  3  dB  bandwidths 
were  2  to  640  Hz  for  the  hydrophone  channels  and  2  to  350  Hz  for  the 
geophone  channels.  Calibration  of  the  system  was  performed  by  command¬ 
ing  a  built-in  autocalibrate  signal. 

The  hard-wired  system  was  very  similar  to  the  buoy  system  except 
the  RF  link  was  replaced  with  a  multi-conductor  cable.  Preamps  were 
mounted  with  the  sensor  as  the  in-water  component  of  the  system. 
Variable-gain  amplifiers  (VGA)  were  mounted  on  the  fixed  oceanographic 
platform  where  the  data  was  recorded.  The  geophone  system  frequency 
response  was  the  same,  while  the  hydrophone  response  was  modified  to 
cover  a  range  of  20  to  4000  Hz.  Commands  were  sent  directly  via  the 
cable  to  change  the  VGA  and  preamp  gain. 
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3.3.2  Sensors .  Two  types  of  sensors  were  connected  to  the 
submerged  unit,  depending  on  the  test  to  be  performed.  One  sensor  was 
the  vertical  array,  while  the  other  was  the  bottom-sensor  unit  consist¬ 
ing  of  a  top-mounted  and  a  bottom-mounted  hydrophone  and  a  three- axis 
geophone.  The  two  systems  were  never  operated  together,  since  a  change 
in  the  internal  electronics  accompanied  the  sensor  change. 

The  vertical  array  used  in  the  tests  off  Panama  City  and  Daytona 
Beach  consisted  of  eight  omnidirectional  LC-10  hydrophones  spaced 
3.7  m  apart.  The  sensitivity  of  each  hydrophone  was  approximately 
-210  dBv/yPa.  In  water  depths  of  31  m,  the  complete  array  was  used, 
and  the  top  of  the  supporting  float  was  approximately  2  m  below  the 
surface.  In  a  water  depth  of  18  m  off  Panama  City,  the  top  three 
hydrophones  were  tied  together  at  the  float  2  m  below  the  surface.  In 
a  water  depth  of  18  m  off  Jacksonville,  the  top  three  hydrophones  were 
tied  back  down  the  array  and  secured  at  the  midway  points  between  the 
lower  hydrophones  as  shown  in  Figure  3. 

In  the  Corpus  Christi  and  Scotian  Shelf  tests,  a  new  array  was 
used  that  contained  removable  omnidirectional  hydrophones  having  the 
same  spacing  and  placement  in  the  water  column  as  the  LC-10  array.  The 
sensitivity  of  each  hydrophone  was  -199  dBv/yPa.  The  upper  hydrophones 
were  tied  back  down  the  array  for  the  shallow-depth  tests  using  this 
array. 

The  hydrophones  used  in  the  bottom-sensor  unit  were  the  top 
hydrophone  with  a  sensitivity  of  -185  dBv/yPa,  and  the  bottom 
hydrophone,  with  a  sensitivity  of  -195  dBv/uPa.  The  top  hydrophone 
was  omnidirectional,  while  the  bottom  hydrophone  was  more  directional 
due  to  its  inset  location  in  the  bottom  bulkhead. 
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The  three-axis  geophone  on  the  bottom-sensor  package  consisted 
of  two  horizontal  and  one  vertical  geophone  mounted  in  three  faces  of 
an  aluminum  cube  bolted  to  the  inside  face  of  the  bottom  bulkhead.  The 
geophones  were  sensitive  to  particle  velocity  in  the  horizontal  or 
vertical  plane.  The  sensitivity  of  the  geophones  at  100  Hz  with  5  dB 
of  preamp  gain  was  2.75  V/cm/sec. 

3.3.3  Sources.  Several  types  of  sources  have  been  used, but  the 
most  widely  used  were  the  7.1  g  (1/4  oz) ,  the  454  g  (1  lb),  and  the 
816  g  [1.8  lb  air-dropped  SUS  (Signal,  Underwater  Sound) ]  charges. 

The  charges  were  composed  of  different  explosive  types,  e.g.,  the  7  g 
charges  were  PETN;  the  weights  given  are  the  equivalent  of  the  same 
weight  of  TNT.  For  all  but  the  largest  charge,  the  detonation  depth 
was  controlled  by  using  a  weighted  charge  at  the  end  of  a  long  cord  of 
the  desired  length,  suspended  from  a  surface  float.  The  SUS  charge  was 
preset  to  detonate  at  a  particular  pressure;  detonation  settings  of  60, 
800,  and  1500  ft  were  available.  Ship-launched  charges  were  rigged 
with  a  long  fuse  to  allow  a  separation  in  range  between  the  ship  and 
the  detonation.  The  physics  of  the  explosive  process  and  the 
characteristics  of  the  signals  were  discussed  by  Weston  (1960). 

In  determining  transmission  loss  in  one-third-octave  bands,  the 
source  level  in  the  same  frequency  bands  was  required.  Since  it  was 
not  possible  to  measure  this  level  for  all  combinations  of  charge 
weight,  depth  and  frequency,  procedures  have  been  devised  by  Gaspin  and 
Shuler  (1971)  for  predicting  these  levels.  Other  researchers  have 
noted  errors  at  low  frequency  in  the  predicted  source  levels  and  have 
published  procedures  for  making  a  low-frequency  correction,  e.g., 

Hughes  (1976). 
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3.4  Analysis  System  and  Procedures 

In  this  section,  the  Hewlett-Packard  Fourier  Analyzer  System  is 
described.  Its  operation  is  related  to  the  calculation  of  the  received 
energy  for  both  a  single  hydrophone  and  an  array  of  hydrophones,  the 
calculation  of  dispersion  within  a  signal,  and  seismic-refraction 
analysis. 

3.4.1  Hewlett-Packard  Fourier  System.  The  HP  system  consists 
of  an  HP  2100  A  computer,  a  hard-wired  Fast  Fourier  Transform  Unit,  an 
oscilloscope  display,  video  terminal,  digital  tape  unit,  disc  memory, 
and  a  control  keyboard.  In  general,  the  system  is  designed  to  perform 
time-and-f requency-domain  analysis  of  data  containing  frequencies  from 
dc  to  100  kHz.  Input  to  the  system  is  either  via  a  two-channel  analog 
input  into  a  A/D  converter  or  from  the  digital  tape.  Upon  keyboard 
command,  the  unit  can  perform  such  operations  as  calculate  forward  and 
inverse  Fourier  transforms,  power  and  cross-power  spectra,  coherence 
functions,  auto-  and  cross-correlation,  ensemble  averaging  in  time  and 
frequency,  etc.  Sampling  and  data-block  size  parameters  are  selectable; 
for  most  of  the  operations  reported  here,  a  sampling  period  of  500  usee 
was  used  which  yields  a  frequency  resolution  of  about  0.5  Hz  for  a 
block  size  of  4096  data  words.  These  parameters  give  a  maximum 
frequency  of  1000  Hz  and  a  maximum  time  window  of  about  2  sec.  All 
input  data  are  calibrated  and  carry  a  scale  factor  describing 
amplitude  and  a  frequency  code  describing  sampling  parameters. 

Keyboard  programs  may  be  written  for  any  sequence  of  steps  involved  in 
a  computation.  BASIC  and  FORTRAN  programs  may  also  be  written  for  the 
HP  2100  A  computer  and  called  as  a  sub-program  from  a  keyboard  program. 
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3.4.2  Received  energy  calculations.  Received  signals  are 
recorded  at  sea  on  analog  tape  along  with  a  time  code.  By  selecting 
the  desired  time  on  a  synchronization  unit,  any  point  on  the  tape  may 
be  selected  within  one  millisecond  in  time.  The  signal  is  low-pass 
filtered  (0  to  1000  Hz)  to  prevent  aliasing  in  the  sampling  operation, 
and  at  the  selected  time,  the  signal  is  strobed  into  the  HP  temporary 
storage  as  a  digital  signal.  The  signal  is  now  ready  for  processing. 

The  purpose  of  processing  the  experimental  data  was  to  obtain 
the  attenuation  of  the  first  mode.  Where  available,  long-range  data 
were  used  for  this  purpose  since,  in  most  cases,  the  energy  was  pre¬ 
dominantly  in  the  first  mode.  Short-range  data  were  also  used,  but 
frequently  these  data  were  corrupted  by  strong  contributions  from  the 
higher-order  modes.  When  it  was  verified  that  only  the  first  mode  was 
present,  the  received  energy  was  calculated  from  a  single  hydrophone 
output.  When  higher-order  modes  were  present,  the  data  were  processed 
using  array  techniques  involving  several  hydrophones . 

The  operations  involved  in  the  single  hydrophone  energy 
calculation  are  diagrammed  in  Figure  4.  The  signal  was  first  filtered 
(anti-aliasing)  and  then  stored  in  the  temporary  storage  of  the  HP. 

The  signal  was  then  transformed  into  the  frequency  domain  by  the 
Fourier  transform  operation  (F) .  The  desired  analysis  filter  was 
obtained  by  calling  for  this  program  as  a  subroutine,  and  the 
frequency-domain  signal  was  then  multiplied  by  the  filter  to  give 
Y(f)  ,  the  filtered  signal.  The  energy  in  the  filtered  signal  was 
obtained  using  the  equation: 
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Figure  4.  Energy  Calculation — Single  Hydrophone 
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Integration  was  performed  as  a  summation  at  each  frequency  (or  channel 
of  the  data  word)  from  f^  to  f^  •  The  energy  levels  were  then 
converted  to  the  logarithmic  domain,  and  the  final  value  in  the  block 
was  printed  to  yield  the  energy  level  in  dB. 

Basically,  two  different  types  of  analysis  "filters"  were  used 
in  processing  the  received  energy.  The  first  was  a  "filter"  approx¬ 
imating  the  well  known  one-third-octave  filter.  As  implemented  on  the 
HP  analyzer,  this  filter  was  approximated  by  summing  the  energy  in 
groups  of  channels  within  the  4096  channel  data  block.  The  frequency 
range  of  each  group  of  channels  represented  approximately  one-third  of 
an  octave  in  frequency.  Nominal  center  frequencies  of  25,  31.5,  40, 

50,  63,  80,  100,  125,  160,  200,  250,  315,  400,  500,  630,  800,  and  1000 
Hz  were  used. 

The  other  type  of  "filter"  is  an  approximation  to  a  finite- 
impulse-response  filter  having  a  Butterworth  amplitude  response  curve. 

In  many  applications,  the  frequency-dependent  phase  response  (or 
envelope  delay)  requires  corrections  when  time  delays  or  velocities  are 
important.  Because  of  this,  filters  have  been  constructed  having 
Butterworth  amplitude  characteristics,  but  with  zero  phase  terms.  This 
is  accomplished  by  setting  the  real  part  of  each  channel  in  the  data 
block  equal  to  the  desired  amplitude  response,  at  that  frequency,  and 
setting  the  imaginary  part  equal  to  zero.  With  this  type  of  filter,  no 
delay  correction  is  necessary.  Low-pass,  three-pole  filters  have  been 
implemented,  the  frequency  and  time  response  tested,  and  the  filters 
used  to  process  dispersion  data.  For  the  array-processed  data  discussed 
next,  bandpass  filters  were  used  having  approximately  one-half-octave 
bandwidths  and  center  frequencies  of  50,  100,  200,  400,  600,  and  750  Hz. 
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As  noted  in  the  previous  paragraph,  the  short-range  data  may 
have  higher-order  modes  present.  One  way  to  extract  the  first-mode 
energy  is  to  use  modified  array  processing,  which  uses  selected 
hydrophones  from  the  array  and  the  modal  distribution  of  pressure  in 
the  water  column  to  obtain  a  spatial  separation  of  the  first  mode. 

This  concept  is  illustrated  in  Figure  5  which  shows  a  typical  modal- 
pressure  distribution  in  the  water  column  for  the  first  two  modes  and 
the  array  hydrophone  locations.  To  cancel  the  second  mode,  hydrophones 
2  and  7  are  added  and,  since  the  phase  of  the  second  mode  differs  by 
180°  for  these  two  signals,  cancellation  will  occur.  An  implementation 
of  this  type  of  processing  where  all  opposite  phase  hydrophone  outputs 
are  summed  is  described  by  Caswell  (1979). 

In  this  paper,  the  two  hydrophones  were  processed  as  indicated 
in  Figure  6.  Signals  A  and  B  were  cross-correlated  to  determine  any 
phase  shift  due  to  tilting  of  the  array  away  from  vertical  in  the  water. 
This  phase  error  was  determined  as  the  number  of  shifts  required  for 
perfect  alignment  (the  maximum  cross-correlation  value) .  Once  this 
was  determined,  each  signal  was  cleared  of  all  the  unwanted  noise. 

One  signal  was  shifted  with  respect  to  the  other  to  correct  the  phase 
error,  and  the  signals  were  then  averaged.  The  remaining  operations 
transformed  the  averaged  signal  back  to  the  frequency  domain  to 
calculate  the  energy. 

3.4.3  Dispersion  analysis.  The  purpose  of  performing  dispersion 
analysis  is  to  infer  characteristics  of  the  ocean  bottom  by  comparing 
experimental  values  of  group  velocity  with  those  predicted  by  3 
normal-mode  model.  Dispersion  is  defined  as  a  variation  in  the 
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Figure  5.  Mode  Amplitude  and  Array  Processing 
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propagation  velocity  with  frequency  and  as  it  results  in  this  case, 
from  the  geometrical  constraints  imposed  by  the  shallow-water  sound 
channel,  it  is  termed  geometrical  dispersion.  Another  tyj e  of 
dispersion,  i.e.,  intrinsic  dispersion,  may  result  for  example,  from 
density  variation  within  a  material.  For  the  shallow-water  sound 
channel,  dispersion  is  controlled  by  the  water  depth,  the  compressional 
sound  speed  in  the  water  and  sub-bottom  layers,  the  layer  thickness  of 
sub-bottom  layers,  and  the  density  in  the  water  and  bottom.  This  type 
of  analysis  was  pioneered  by  Pekeris  (1948)  using  the  experimental 
results  of  Worzel  and  Ewing  (1948) . 

Dispersive  propagation  for  a  water  layer  and  a  very  thick 
bottom  layer  is  illustrated  in  Figure  7  for  the  first  mode.  Part  (a) 
shows  a  typical  plot  of  the  group  velocity  as  a  function  of  frequency, 
and  part  (b)  illustrates  how  this  behavior  occurs  in  a  typical, 
although  idealized,  shallow-water  received  signal.  As  can  be  seen,  the 
signal  starts  with  a  low-frequency  (f^)  refracted  arrival.  At  a  time 
corresponding  to  the  sound-speed  in  water  (c) ,  a  high-frequency 
watet-wave  arrival  is  superimposed  upon  the  low-frequency  rider  wave 
(f  ).  The  frequency  of  the  water  wave  gradually  decreases  until  the 
Airy-phase  frequency  (f^)  is  reached.  The  frequency  then  remains 
constant  until  the  signal  decays  into  the  background  noise.  The 
refracted  low-frequency  arrival  travels  with  the  velocity  of  the 
bottom  layer  and  has  a  characteristic  frequency  termed  the  cutoff 
frequency. 

The  experimental  dispersion  data  was  obtained  by  first  filtering 
the  signal  with  several  zero-phase-shift,  low-pass  filters  having 
different  cutoff  frequencies  (e.g.,  25,  50,  100,  200  Hz).  Each 


Figure  7.  Geometrical  Dispersion  in  Shallow  Water 
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filtered  signal  was  plotted  in  the  time  domain  and  the  arrival  times 

(t  )  of  signal  peaks  were  listed.  These  were  measured  with  respect  to 
n 

the  water-wave  arrival  time  since  this  velocity  was  known.  The 
frequency  of  the  nC^  arrival  was  calculated  from  the  relation 


the  group  velocity  was  calculated  from  the  equation 


n  R  +  t  c 
—  n 

The  times,  t  ,  were  measured  with  respect  to  the  water-wave  arrival, 
n 

for  the  range  R  and  the  known  sound-speed  in  the  water  c  .  The 
plus  (+)  sign  was  used  for  arrivals  after  the  water  wave,  and  the  minus 
(-)  sign  was  used  for  the  refracted  arrivals  occurring  before  the  water 
wave.  This  procedure  was  followed  for  signals  from  several  ranges  and 
for  each  analysis  filter  used.  If  the  water  and  bottom  characteristics 
were  independent  of  range,  the  plotted  data  for  shots  from  different 
ranges  should  fall  along  the  same  curve. 

Once  the  experimental  values  of  f^  and  u^  were  determined, 
the  analysis  was  continued  by  assuming  values  for  the  compressional 
sound  speeds,  layer  thicknesses,  and  densities  in  the  sub-bottom  layer; 
it  was  assumed  that  the  water  velocity  and  depth  were  known.  These 
parameters  were  then  used  as  inputs  to  a  normal-mode  propagation  model 
from  which  predicted  values  of  group  velocity  were  obtained  as  a 
function  of  frequency.  These  predicted  values  were  then  compared  with 
the  experimental  values,  the  assumed  bottom  parameters  changed,  the 
model  run  again,  and  new  values  for  the  group  velocity  predicted  for 
another  comparison.  The  procedure  was  repeated  until  the  predicted 


53 


values  matched  the  experimental  values  with  reasonable  accuracy. 

Once  this  occurred,  the  most  recent  bottom  model  used  was  assumed 
to  be  the  correct  bottom  model. 

3.4.4  Seismic-refraction  analysis.  The  use  of  dispersion 
analysis  for  determining  ocean-bottom  parameters  has  been  discussed. 
However,  an  initial  trial  model  must  be  provided  as  a  starting  point 
for  the  dispersion  analysis.  Frequently,  this  initial  model  is  based 
on  seismic-refraction  analysis.  This  type  of  analysis  relies  on  the 
fact  that  when  sound  rays  strike  the  ocean  bottom  at  the  critical 
angle  (assuming  all  velocities  are  increasing  with  depth) ,  energy 
travels  along  the  ocean-bottom  interface  with  a  velocity  equal  to  the 
compressional  sound  speed  in  the  bottom.  Energy  also  penetrates  to 
deeper  layers,  so  that  critical-angle  refraction  occurs  at  each 
velocity  discontinuity  (i.e.,  layer)  found  in  the  bottom.  These 
interface  waves  are  eventually  refracted  back  toward  the  water  layer 
where  they  may  be  received  on  a  hydrophone.  As  noted,  the  key 
assumption  is  that  the  velocity  in  successively  deeper  layers  is 
always  increasing. 

For  the  seismic-refraction  experiments,  charges  were  detonated 
on  the  bottom  at  0.23  km  range  intervals  for  the  first  1.9  km  and 
then  every  0.46  km  out  to  the  maximum  range — usually  8  to  12  km. 

Long  fuses  were  employed  to  insure  bottom  detonation.  This  required 
a  range  correction  for  the  fuse-burn  time.  Ranges  were  calculated 
based  on  the  signal  received  on  a  hydrophone  deployed  from  the  shot 
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The  theory  and  analysis  of  the  refracted  arrivals  can  be 

understood  with  the  aid  of  the  four-layer  model  in  Figure  8.  The 

model  represents  shots  detonated  on  the  seabed  at  point  and 

received  on  the  seabed  at  point  .  At  the  receiver,  assuming  a 

large  separation  R  ,  the  first  signal  arrives  after  traveling  with 

time  t^  >  the  next  travels  for  t^  seconds,  then  t^  ,  and  finally, 

at  t  ,  the  water  wave  arrives.  Implied  in  this  sequence  is  that 

c.  >  c„  >  c,  >  c  as  stated  earlier.  With  the  exception  of  t  ,  all 
3  2  1  w  r  w 

are  refracted  arrivals  that  obey  the  travel-time  relations  in  the  lower 
part  of  the  figure.  Each  of  these  is  derived  using  Snell's  law  and 
geometrical  ray  theory.  Snell's  law  is  stated  for  this  situation  by 
the  two  assumptions  in  the  figure,  where  it  is  implied  that  the 
refracted  ray  leaves  the  interface  at  an  angle  of  90°;  that  is, 

Sln  all  =  sin  90°  =  _1_ 

C1  C2  C2 

Note  that  t^  ,  c2  ’  an<*  C3  are  aH  linear  functions  of  R  .  Thus, 
t^  as  a  function  of  R  is  a  straight  line  with  a  slope  of  1/c^  and 
a  zero  intercept.  A  plot  of  t ^  as  a  function  of  R  has  a  slope  of 
I/C2  and  a  time-axis  intercept  of  (2h^  *c2  ~  cl~^clc2  ’  C3  °^e>’s 

a  similar  relationship.  This  is  illustrated  in  Figure  9. 

With  this  background,  the  analysis  procedure  can  be  described 
as  follows.  The  received  signals  were  recorded  with  a  high  gain  so 
that  the  refracted  arrivals  were  clearly  visible.  Upon  playback  using 
the  HP,  they  were  plotted  with  an  expanded  scale  so  that  the  time  of 
arrival  could  be  measured  with  some  accuracy.  The  first-arrival  points 
were  recorded  and  plotted  for  each  shot.  Since  the  range  increased  for 
each  shot,  these  points  were  plotted  as  the  crosses  in  Figure  9.  Once 
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the  data  was  all  plotted,  points  lying  on  the  same  line  were  obvious 
and  could  then  be  grouped  together.  A  straight-line-fitting  procedure 
was  applied  to  each  of  these  groups  to  obtain  the  slopes  and  intercepts. 
Thus,  from  the  line  identified  as  t^  ,  the  velocity  c^  was  obtained, 
and  knowing  the  intercept  value,  h^  was  readily  calculated.  Finally, 
from  t^  ,  the  velocity  c^  and  thickness  were  calculated.  This 

procedure  may  be  extended  to  any  number  of  layers  as  noted  by  Officer 
(1958). 

In  reality,  certain  bottom  structures  require  a  more  detailed 
analysis  or  may  give  incorrect  results.  For  example,  a  thin  layer  may 
not  appear  as  a  first  arrival,  so  frequently  the  arrival  times  of  the 
secondary  arrivals  must  be  obtained  as  well — in  most  cases,  a  difficult 
task.  If  a  low-velocity  layer  (a  layer  with  a  velocity  less  than  the 
layer  above  it)  occurs,  all  layer  thicknesses  calculated  below  the 
low-velocity  layer  will  be  in  error.  This  problem  is  discussed  by 
Dobrin  (1976).  We  have  assumed  a  flat,  horizontal  bottom  structure; 
however,  sloping  bottoms  may  be  analyzed  using  these  techniques  as  well. 
In  this  case,  the  inverse  slope  is  no  longer  the  velocity,  but  is  a 
function  of  the  critical  angle  and  the  angle  of  the  bottom  slope.  A 
modification  to  the  experimental  geometry,  i.e.,  switching  the  source 
and  receiver  position,  is  frequently  employed  to  resolve  the  unknown 
slope  angle.  Both  Officer  and  Dobrin  discuss  this  problem. 

In  general,  seismic-refraction  analysis  gives  good  results 
where  velocities  are  concerned,  but  has  limits  on  the  order  of  a 
wavelength  or  so  for  thickness  measurements.  Dispersion  analysis  is 
employed  to  give  better  resolution  for  detecting  thin  near-surface 
bottom  layers. 
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3.5  Data  Characterizing;  the  Sediments 

Several  different  types  of  input  data  are  required  in  order  to 
characterize  the  bottom  for  sediment-modeling  purposes.  The  dispersion 
and  seismic-refraction  analyses  are  used  to  describe  a  large  section  of 
the  seabed.  The  sediment  model  (discussed  in  Chapter  II)  may  describe 
only  the  top  few  meters  or  so.  The  data  required  for  this  is  grain- 
size  and  distribution,  bulk  density,  porosity,  etc.,  which  may  be 
obtained  from  surface  grab-samples  or  cores.  In  certain  cases, 
specialized  acoustic-profiling  systems  (e.g.,  the  Huntec  DTS  System) 
can  give  some  measure  of  the  grain-size  and  distribution.  With  this 
data,  as  a  minimum,  the  mean  grain-size  and  sorting  (i.e.  ,  standard 
deviation) ,  the  permeability,  and  pore-size  parameter  may  be  calculated 
for  use  in  the  model.  The  bulk  density  may  be  calculated  once  the 
porosity  is  estimated  from  empirical  data.  The  desired  output  of  the 
model  is  the  sediment  absorption  as  a  function  of  frequency. 

In  all  cases,  attempts  were  made  to  obtain  core  and  grab-sample 
data  as  part  of  the  experimental  program.  In  most  cases,  however,  the 
available  equipment  was  inadequate  for  obtaining  data  in  the  bottom 
types  encountered.  Because  of  this,  core  and  grab-sample  data 
collected  by  other  researchers  were  frequently  obtained  by  contractual 
analysis.  Finally,  where  no  other  data  could  be  obtained,  the  open 
literature  was  studied  for  the  best  available  data.  The  availability, 
quantity,  and  type  of  sediment  data  are  discussed  on  a  site-by-site 
basis  in  the  appropriate  chapters. 


CHAPTER  IV 


BOTTOM- ABSORPTION  ANALYSIS 
DAYTONA  BEACH  EXPERIMENTS 


4.1  Introduction 

Experimental  data  obtained  at  two  sites  off  the  coast  of  Daytona 
Beach,  Florida,  is  analyzed  to  determine  the  magnitude  of  losses 
attributable  to  absorption  in  the  sedimentary  layer.  The  experiments 
performed  are  discussed  in  detail  in  Chapter  III,  but  are  briefly 
reviewed  here  to  establish  the  applicable  environmental  and  experimental 
parameters  for  the  test  sites  under  discussion.  The  analysis  procedures 
discussed  in  Chapter  III  are  applied  to  recorded  data  to  obtain  mode- 
attenuation  coefficients  from  which  the  experimental  bottom-absorption 
losses  are  determined.  The  bottom  models  established  for  both  sites 
and  the  techniques  used  to  obtain  them  are  discussed.  Finally,  the 
sediment  data  available  for  each  site  is  presented,  and  the  bottom 
absorption  predicted  by  the  sediment  models  is  compared  with  the 
experimental  values. 

4. 2  Experiments 

Experiments  were  performed  during  November  of  1976  at  two  sites, 
depicted  in  Figure  10,  off  the  coast  of  Daytona  Beach,  Florida.  The 
water  depth  was  nominally  19  m  at  the  shallow  site  (Site  1)  and  32  m 
at  the  deeper  site  (Site  2)  .  In  addition  to  the  data  analyzed  from  the 
1976  experiment,  certain  data  obtained  in  the  vicinity  of  Site  2  during 
an  earlier  experiment  (1973)  were  analyzed  as  well.  Sound-velocity  and 


Figure  10.  Daytona  Beach  Test  Area 
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XBT  data  were  obtained  at  the  sites  and  along  the  propagation  tracks. 

A  12  kHz  fathometer  was  also  operated  along  the  tracks  to  obtain  the 
bathymetry.  Neither  cores  nor  sub-bottom  profiling  data  were  obtained 
during  the  tests.  Because  of  the  lack  of  bottom  data,  grab-samples  and 
cores  were  obtained  from  sources  to  be  discussed  later. 

At  Site  1,  the  experimental  data  obtained  were  from  the  mode- 
attenuation  and  seismic-refraction  experiments.  Sources  for  the 
mode-attenuation  tests  were  7.1  g  charges  detonated  at  depths  of  9.5 
and  12.8  m.  The  receiver  was  the  vertical  array  and  since  the  array 
length  was  greater  than  the  water  depth,  it  was  folded  back  (see 
Figure  3)  to  provide  closer  spacing  in  the  upper  half  of  the  water 
column.  The  variation  in  water  depth  and  sound  velocity  with  range 
along  the  constant-depth  (Leg  E)  and  deep-to-shallow  (Leg  F)  legs  are 
depicted  in  Figure  11.  The  long-range  transmission  loss  experiment  was 
not  performed  at  this  site  due  to  adverse  weather  conditions.  Signals 
selected  from  the  constant-depth  ^eg  shots  of  the  seismic-refraction 
experiment  were  used  for  dispersion  analysis  at  this  site. 

The  experiments  performed  at  Site  2  were  very  similar  to  those 
performed  at  Site  1.  Short-range,  mode-attenuation  data  were  obtained 
at  Site  2  along  both  sloping  (Leg  H)  and  constant-depth  (Leg  G)  tracks 
out  to  a  range  of  12  km.  The  water  depth  shown  in  Figure  12  was 
essentially  constant  over  both  legs.  The  small  7.1  g  charges  were 
detonated  at  8.5  and  17  m  and  received  on  the  same  vertical  array 
extended  to  its  full  length.  Sound-velocity  profiles  taken  indicated 
that  the  water  was  isovelocity. 

Long-range  propagation  data  were  obtained  at  Site  2  out  to  a 
range  of  95  km  on  the  constant-depth  (Leg  J)  track  depicted  in 
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Figure  11.  Depth  and  Sound  Velocity  Profiles — Daytona 
Beach  Site  1 


Figure  12.  Air-dropped  bathythermograph  probes  (AXBT's)  indicated 
that  temperatures  varied  little  with  water  depth. 
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Seisnic-refraction  experiments  were  performed  at  Site  2,  during 
both  the  1976  and  the  1973  tests  .  Also  performed  during  the  1973 
tests  was  a  short-range  seismic  experiment  using  an  array  of  12 
hydrophones  placed  in  a  straight  line  in  contact  with  the  bottom. 

Small  7.1  g  charges  were  detonated  on  the  bottom  at  approximately  31  m 
from  the  closest  hydrophone,  and  arrival  times  were  measured  to  calculate 
refracted-wave  velocities.  The  water  depth  for  this  experiment  was 
43  m. 

4.3  Analysis  of  Data — Site  1 

The  objective  of  the  data  analysis  for  both  Sites  1  and  2  was  to 
show  that  bottom  absorption  can  be  predicted  by  the  sediment  model.  To 
show  this,  experimental  values  of  bottom  absorption  were  calculated  as 
a  basis  for  comparison.  This  was  accomplished  by  first  calculating  the 
received  energy  and  then  the  mode-attenuation  coefficients  for  the 
first  mode  from  which  the  experimental  bottom-absorption  coefficients 
were  obtained. 

Since  long-range  data  were  not  available  for  Site  1,  short-range 
data  were  used.  Hance,  some  consideration  for  higher-order  modes  was 
required  and,  consequently,  modified  array  processing  of  the  received 
signals  (described  in  Chapter  III)  was  performed.  Signals  received 
from  the  small  charges  detonated  at  9.5  m  were  analyzed  by  these 
techniques . 

4.3.1  Received  energy  and  mode-attenuation  calculations — Site  1. 
Short-range  data  were  obtained  along  Leg  E  for  ranges  of  0.76  to  12.0 


t 
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km  and  along  Leg  F  for  0.5  to  15.9  km.  Two  shots  were  detonated  at 
the  9.5  m  source  depth  at  each  range  point.  Data  were  obtained  over 
a  frequency  range  of  50  to  750  Hz. 

To  obtain  the  received  energy  due  only  to  the  first  mode,  a 
summation  of  two  hydrophones  was  performed  in  the  laboratory  in  order 
to  cancel  the  contribution  of  the  higher  modes.  This  calculation  was 
performed  by  first  selecting  the  two  hydrophone  signals  having  the  best 
signal-to-noise  ratio  and  location  in  the  water  column  to  obtain 
optimum  cancellation  of  the  second  mode.  For  Leg  E,  the  two  hydrophones 
used  were  numbers  2  and  6,  and  for  Leg  F,  numbers  2  and  7  were  used. 
Analysis  filters  (one-half  octave)  centered  at  50,  100,  200,  400,  600, 
and  750  Hz  were  used.  The  hydrophone  pairs  were  processed  as  noted  in 
Section  3.4.2  to  obtain  the  received  energy  for  the  first  mode. 

In  order  to  calculate  the  mode-attenuation  coefficients,  the 

received  energy  was  corrected  for  system  gain  and  cylindrical  spreading. 

Since  the  mode-attenuation  coefficients  were  defined  as  the  slope  of 

the  line  obtained  by  plotting  the  loss  as  a  function  of  range,  only 

range-dependent  effects  were  considered.  Hence,  the  source  strength 

was  not  required  for  this  calculation  and  the  transmission  loss  was 

not  calculated.  The  received  energy  was  corrected  for  gain  relative 

to  the  gain  at  one  range  and  corrected  for  cylindrical  spreading 

relative  to  the  spreading  loss  at  one  range.  Data  were  deleted  for 

those  shots  where  the  signal-to-ncise  ratio  was  less  than  approximately 

3  dB.  The  zero  decibel  reference  for  the  relative  energy  calculation 

2 

is  defined  as  1  volt  -sec. 

Using  the  relative  energy  levels,  the  mode-attenuation 
coefficients  were  obtained  by  performing  a  straight-line,  least  square 
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fit  to  the  data  at  each  frequency.  The  mode-attenuation  coefficients 
are  shown  in  Table  3  for  Legs  E  and  F  along  with  the  95  percent 
confidence  interval  limits.  The  relative  energy  levels  and  straight 
line  fits  are  plotted  in  Figures  13  and  14  for  Legs  E  and  F, 
respectively. 

4.3.2  Bottom  model — Site  1.  To  calculate  the  bottom 
absorption  given  the  experimental  mode-attenuation  coefficient  for  the 
first  mode  (6),  the  energy  present  in  the  sub-bottom  layers  (y)  must 
be  known.  The  normal-mode  program  gives  this  information;  however,  to 
run  the  program,  a  valid  model  for  the  sub-bottom  structure  must  be 
obtained  for  use  as  input  to  the  program. 

Such  a  sub-bottom  model  has  been  determined  by  an  analysis  of 
the  dispersion  exhibited  by  shots  received  along  the  constant-depth 
track  (as  discussed  in  Chapter  III).  The  experimental  dispersion  data 
is  plotted  in  Figure  15  as  the  discrete  points.  Group  velocities 
predicted  by  a  four-bottom-layer  normal-mode  propagation  model,  using 
the  water-bottom  structure  shown,  are  plotted  as  the  solid  line. 

Although  the  water  depth  varied  along  Legs  E  and  F,  as  shown  in 
Figure  11,  a  water  depth  of  20.5  m  was  used  for  both  legs.  The 
sound-velocity  profile  shown  is  from  a  range  of  3.8  kra  along  Leg  F 
and  is  representative  of  those  encountered  along  the  track.  The  use 
of  this  profile  for  Leg  E  as  well  has  little  effect  on  the  comparison 
of  predicted  and  experimental  data. 

The  bottom  model  shown  in  Figure  15  was  selected  on  the  basis  of 
dispersion  analysis  and  is  consistent  with  the  results  of  other 
researchers  who  have  studied  the  test  area.  Meisburger  and  Field  (1975) 
have  conducted  detailed  surveys  over  the  inner  continental  shelf  from 
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TABLE  3 


MODE- ATTENUATION  COEFFICIENT  (6)— SITE  1 
Daytona  Beach 


Mode-Attenuation 

95%  Confidence 

Interval 

Frequency 

(Hz) 

Coefficient  (6) 
(dB/km) 

Upper  Limit 

Lower  Limit 

50 

3.09 

3.99 

2.19 

100 

1.29 

1.59 

0.99 

200 

1.33 

1.75 

0.91 

400 

1.08 

1.39 

0.77 

600 

1.07 

1.60 

0.54 

750 

1.43 

1.71 

1.15 

50 

2.78 

3.35 

2.21 

100 

1.46 

1.86 

1.06 

200 

1.27 

1.53 

1.01 

400 

0.97 

1.31 

0.63 

600 

1.03 

1.41 

0.  o5 

750 


1.38 


1.66 


1.10 


6 


8 

RANGE  ( km ) 


ilative  Energy  Levels — 
2ach  Sice  1,  Leg  F 


GROUP  VELOCITY  (m /sec) 


Figure  15.  Daytona  Beach  Experimental  and 
Model  Dispersion — Site  1 
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Cape  Canaveral  to  Georgia.  Many  grab-samples  and  cores  were  obtained 
to  correlate  with  seismic-reflection  lines  run  perpendicular  to  the 
shore.  Two  of  these  reflection  lines  bound  Site  1  to  the  north  and 
south  and  show  primary  reflectors  (i.e.,  a  layered  structure)  consistent 
with  the  bottom  model  of  Figure  15.  The  velocities  used  in  the  present 
bottom  model  show  reasonable  agreement  with  seismic-refraction  results 
obtained  by  Kersey  et  al.  (1950),  Antoine  and  Henry  (1965),  and 
Sheridan  et  al.  (1966)  off  the  Florida  coast.  The  uppermost  layer  is 
reported  to  have  a  velocity  between  1630  and  1810  m/sec  for  all  of  the 
results  reported;  the  next  layer  has  a  velocity  range  of  2240  to 
2690  m/sec.  The  layer  thickness  obtained  by  seismic-refraction 
analysis  are  generally  much  thicker  than  found  by  Meisburger  and  Field, 
and  probably  reflect  the  general  insensitivity  of  the  seismic-refraction 
techniques  for  showing  small-scale  changes  in  thickness. 

4.3.3  Experimental  sediment  absorption — Site  1.  With  the 

four-layer  normal-mode  program  and  four-layer  bottom  model  of  Figure  15, 

mode-amplitude  distributions  of  pressure  throughout  the  water  column 

and  bottom  were  calculated  and  values  for  y  in  each  sub-bottom  layer 

obtained.  This  data  is  listed  in  Table  4  where  y^  is  the  energy 

(2) 

ratio  for  the  first  bottom  layer,  and  y  is  the  ratio  for  the 
second,  and  so  forth.  The  sum  of  the  three  is  y  .  The  bottom- 
absorption  coefficients  for  Leg  E  and  Leg  F  were  calculated  by 

dividing  the  respective  mode-attenuation  coefficients  by  y  .  These 
values  of  a.  will  be  compared  with  the  bottom  absorption  predicted 


by  the  sediment  model. 
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4.3.4  Predicted  sediment  absorption — Site  1.  As  discussed  in 

Chapter  II,  the  model  for  a  porous  fluid-filled  medium  as  developed  by 

Biot  and  later  refined  by  Stoll  was  used  to  predict  sediment  absorption 

and  velocities.  In  order  to  apply  this  model,  certain  parameters 

describing  the  sediments  in  the  test  area  are  required.  Ideally,  the 

area  should  be  well-sampled  with  deep  cores  analyzed  to  yield  such 

parameters  as  mean  grain-size,  grain-size  statistics,  density,  porosity, 

and  mineralogy.  At  Site  1,  both  cores  and  grab-samples  are  available. 

The  cores  are  numbers  139,  140,  147,  148,  and  173,  discussed  by 

Meisburger  and  Field  (1975) ,  analyzed  at  the  University  of  Texas 

at  Arlington.  The  grab-samples  (65  total)  are  reported  by  Milliman 

et  al.  (1972)  and  were  analyzed  at  Duke  University.  The  grab-samples 

provide  the  best  coverage  near  Site  1  as  may  be  seen  in  Figure  10  where 

both  the  cores  and  grab-samples  are  noted.  The  grab-sample  analysis 

provided  only  the  grain-size,  grain-size  statistics,  and  calcium 

carbonate  content.  However,  empirical  studies,  such  as  performed  by 

Hamilton  (1969) ,  may  be  used  to  relate  grain-size  to  porosity  and  then 

to  density.  The  pertinent  grab-sample  data  is  listed  in  Table  5.  As 

can  be  seen  in  Figure  10,  the  northern  and  southern  sample  lines  bound 

the  site  and  propagation  tracks  on  the  north  and  south.  A  study  of 

the  mean  grain-size  (M  )  of  these  samples  reveals  a  gradual  decrease 

z 

in  grain-size  until  a  minimum  is  reached  at  sample  1482  along  the 
north  and  sample  1478  along  the  south.  The  mean  grain-sizes  increase 
from  the  minima  to  the  end  of  the  sample  lines.  This  behavior  seems  to 
indicate  that  a  band  of  moderately-sorted,  coarser  sediments  runs 
parallel  to  the  shore  across  the  site  and  both  propagation  tracks. 


TABLE  5 
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Note:  phi  =  -  log  d;  J  is  the  mean  grain  diameter  in  mm. 


The  two  sample  lines  are  characterized  for  modeling  purposes  by  the 
average  values  noted  in  Table  5. 

Based  on  these  average  mean  grain-sizes  and  standard  deviations, 

permeability  and  pore-size  parameters  (K  and  a)  were  calculated  as 

noted  in  Chapter  II  for  use  in  the  sediment  model.  The  aggregate  bulk 

modulus  k  for  the  minerals  within  the  sediment  frame  was  calculated 
r 

using  the  VRH  method  (described  in  Chapter  II)  for  a  calcium  carbonate 
content  of  30  percent  obtained  from  Table  5  and  a  70  percent  quartz 
sand  content.  Densities  and  porosities  were  estimated  as  noted  earlier. 
The  parameters  describing  the  two  sediment  models  used  are  listed  in 
Table  6,  and  the  resulting  bottom  absorption  for  the  frequencies  of 
interest  are  listed  in  Table  7. 

4.3.5  Discussion  of  results — Site  1.  The  bottom  absorption 
obtained  from  the  mode-attenuation  coefficients  (the  experimental  data) 
and  that  predicted  by  the  sediment  models  of  Biot  and  Hamilton  are 
compared  in  Figure  16.  For  the  Biot  model,  good  agreement  is  obtained 
between  100  and  400  Hz  for  the  model  representing  a  sediment  with  a 
mean  grain-size  of  1.44  phi.  Much  less  attenuation  is  predicted  by 
the  model  at  50  Hz  than  was  measured.  This  occurs  because  the  sediment 
model  represents  only  a  very  shallow  bottom  layer,  and,  at  50  Hz,  the 
sound  waves  penetrate  into  the  third  bottom  layer.  The  additional 
attenuation  due  to  these  deeper  layers  results  in  the  increased  low- 
frequency  attenuation.  The  Hamilton  model  predicts  greater  absorption 
at  low  frequency  but  less  at  high  frequency,  giving  somewhat  poorer 
agreement  than  the  Biot  model. 


TABLE 


10, 000 


Figure  16.  Comparison  of  Experimental  and 
Predicted  Bottom  Absorption — 
Davtona  Beach  Site  1 
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At  600  and  750  Hz,  the  modal  attenuation  in  Table  3  is  greater 
than  expected,  based  on  the  general  decrease  of  the  coefficient  between 
50  and  400  Hz.  One  possibility  for  this  large  attenuation  coefficient 
is  that  higher-order  modes  are  present  at  the  shorter  ranges  and  a 
gradual  attenuation  (or  stripping)  of  these  modes  occurs  resulting  in  a 
higher  effective  attenuation  for  the  first  mode  as  the  range  increases. 
This  possibility  was  investigated  by  further  processing  on  the 
hydrophone  outputs  which  appeared  to  have  the  least  contribution  due 
to  the  higher  modes.  Typical  750  Hz  signals  are  illustrated  in  Figure 
17  where  the  upper  signal  has  a  great  deal  more  energy  in  modes  two  and 
three  than  does  the  lower  signal.  In  order  to  obtain  the  energy  due 
only  to  the  first  mode,  the  signals  were  processed  by  removing  all  of 
the  signal  except  the  first  mode  due  to  the  shock  pulse.  That  is,  only 
the  signal  between  A  and  B  in  Figure  17  was  retained.  The  result  was 
that  the  mode-attenuation  coefficient  at  750  Hz  was  higher  still.  Thus, 
it  has  been  concluded  that  mode  stripping  is  not  responsible  for  the 
high  mode-attenuation  coefficients  at  600  and  750  Hz. 

Another  loss  mechanism,  associated  with  the  higher  frequencies, 
which  was  considered  as  a  possible  source  of  the  600  to  750  Hz  attenua¬ 
tion  is  scattering  due  to  any  roughness  associated  with  the  surface  and 
bottom.  This  mechanism  was  ruled  out,  however,  since  the  seas  were 
calm  (sea  state  one  or  less)  on  the  particular  day  of  testing  and  the 
bottom  variations  change  slowly  with  range. 

The  last  source  of  attenuation  to  be  examined  as  a  possibility 
was  attenuation  of  a  biological  origin.  This  possibility  was  first 
considered  when  large  masses  of  somewhat  random  reflectors  were  noted 
on  the  12  kHz  fathometer  records  obtained  ns  the  ship  dropped  shots 
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along  Legs  E  and  F.  An  example  of  this  Is  shown  In  Figure  18.  It  Is 
assumed  that  the  reflectors  are  some  type  of  fish  as  they  are  present 
In  varying  degree  along  both  Legs  E  and  F.  This  source  of  attenuation 
has  been  considered  by  Weston  and  Ching  (1970).  They  note  a  diurnal 
pattern  in  the  magnitude  of  attenuation  with  the  maximum  occurring  at 
night.  Attenuation  of  up  to  2  dB/km  at  700  Hz  was  observed  for  a 
concentration  of  one  24-cm  fish  to  approximately  10  square  meters  of 
sea  surface.  If  it  is  assumed  that  the  Biot  sediment  model  accurately 
predicts  the  bottom  absorption  at  750  Hz  [i.e.,  6  =  (0.0009)(707  dB/km 
*  0.64  dB/km],  then  the  excess  attenuation  attributed  to  other 
mechanisms  must  be,  substractlng  this  6  from  that  in  Table  3, 
approximately  0.79  dB/km.  This  is  well  within  the  magnitude  observed 
by  Weston  and  Ching.  However,  whether  the  concentration  and  size  of 
the  scatterers  are  consistent  with  the  observed  attenuation  is  an 
unanswered  question. 

4.4  Analysis  of  Data — Site  2 

At  Site  2,  both  long-range  and  short-range  data  were  obtained. 
The  short-range  data  were  analyzed  by  McDaniel  (1979) .  Mode-attenua¬ 
tion  coefficients  were  calculated  in  that  paper  and  are  presented  for 
both  sloping  (H)  and  constant  depth  (G)  tracks.  Results  from  the 
constant-depth  leg  of  the  long-range  propagation  experiment  are 
presented  in  the  next  section.  A  bottom  model  is  described  based  on 
seismic-refraction  and  dispersion  analysis,  and  the  results  from  the 
short-range  seismic  experiment  conducted  in  1973. 
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Figure  18.  Fathometer  Record — Daytona  Beach  Site  1, 
Leg  E  (5.6  -  7.4  km) 
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4.4.1  Received  energy  and  mode-attenuation  calculations — Site  2. 
An  analysis  of  the  short-range  data  by  McDaniel  indicated  that,  for 
signals  received  on  the  vertical  array  from  the  small  charges  detonated 
at  17  m,  sufficient  suppression  of  the  second  mode  was  obtained.  By 
selecting  the  array  hydrophone  (number  5  from  the  bottom)  nearest  the 
node  for  the  third  mode,  excellent  suppression  of  both  modes  two  and 
three  was  obtained.  The  data  was  then  processed  as  described  in 
Chapter  111  for  the  calculation  of  the  received  energy  using  a  single 
hydrophone.  The  mode-attenuation  coefficients  for  the  first  mode  were 
calculated  by  McDaniel  for  one-third-octave  center  frequencies  of  40  to 
500  Hz  and  are  presented  in  Table  8  along  with  the  95  percent  confidence 
limits  on  the  coefficients. 

The  signals  received  along  the  constant-depth  leg  (Leg  J)  of  the 
long-range  propagation  experiment  were  analyzed  in  one-third-octave 
bands  with  center  frequencies  from  160  to  500  Hz.  The  received  energy 
was  calculated  from  signals  received  on  the  top  hydrophone  of  the  bottom 
unit  (one  meter  above  the  bottom)  according  to  the  procedure  described 
in  Chapter  III  for  a  single  hydrophone.  The  received  energy  was  then 
corrected  for  relative  gain  and  cylindrical  spreading  and  is  plotted  in 
Figure  19.  The  straight  lines  shown  represent  a  least-squares  straight- 
line  fit  to  the  data  points.  The  mode-attenuation  coefficients  (for 
the  first  mode)  obtained  from  the  straight  lines  are  listed  in  Table  9. 

4.4.2  Bottom  model — Site  2.  A  bottom  model  for  Site  2  was 
constructed  from  seismic-refraction  experiments  performed  during  the 
1973  and  1976  tests,  the  short-range  seismic  experiment  performed  in 
1973,  and  dispersion  analysis  performed  on  shots  from  the  1976  tests. 
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TABLE  8 

MODE-ATTENUATION  COEFFICIENT  (6) —SITE  2,  LEGS  G  AND  H 

Daytona  Beach 


Mode-Attenuation 

95%  Confidence 

Interval 

Frequency 

(Hz) 

Coefficient  (6) 
(dB/km) 

Upper  Limit 

Lower  Limit 

40 

2.60 

2.93 

2.27 

63 

1.78 

2.07 

1.49 

100 

1.28 

1.55 

0.98 

160 

0.92 

1.05 

0.79 

200 

0.80 

0.85 

0.75 

250 

0.62 

0.67 

0.58 

315 

0.50 

0.73 

0.28 

400 

0.40 

0.52 

0.28 

500 

0.32 

0.47 

0.18 

40 

5.28 

6.65 

3.91 

63 

4.97 

6.02 

3.92 

100 

2.91 

3.19 

2.64 

160 

1.11 

1.21 

1.00 

200 

0.75 

0.88 

0.63 

250 

0.78 

0.94 

0.62 

315 

0.41 

0.65 

0.16 

400 

0.40 

0.52 

0.28 

500 


0.32 


0.47 


0.18 
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TABLE  9 

MODE -ATTENUATION  COEFFICIENT  (5)— SITE  2,  LEG  J 
Daytona  Beach 

Mode -Attenuation  95%  Confidence  Interval 


Frequency  Coefficient  (6) 


(Hz) 

(dB/km) 

Upper  Limit 

Lover  Limit 

160 

0.75 

0.92 

0.58 

200 

0.65 

0.76 

0.54 

250 

0.50 

0.57 

0.43 

315 

0.43 

0.48 

0.38 

400 

0.38 

0.43 

0.33 

500 

0.34 

0.39 

0.28 

The  seismic-refraction  results  are  summarized  in  Table  10  where  they 
are  also  compared  with  work  by  Hersey  et  al.  (1950),  and  Worzel  and 
Ewing  (1948) .  These  models  were  tested  through  the  use  of  dispersion 
analysis  and  the  results  are  plotted  in  Figure  20.  The  three- layer 
model  shown  in  the  figure  gave  the  best  match  to  the  experimental  data 
and,  although  this  model  differed  somewhat  from  the  other  results  in 
that  the  first  layer  thickness  was  much  less,  it  was  confirmed  by  the 
results  from  the  short-range  seismic  experiment  [described  by  McDaniel 
(1980)].  Such  an  experiment  provided  better  resolution  for  detecting 
thin  layers  than  experiments  using  larger  charges  and  greater  source- 
receiver  separations. 

4.4.3  Experimental  sediment  absorption — Site  2.  Once  a  bottom 
model  is  established,  it  may  be  used  in  the  normal-mode  model  and  the 
distribution  of  energy  among  the  bottom  layers  determined.  This  was 
done  using  the  three-layer  bottom  model  in  Figure  20.  By  assuming  the 
absorption  coefficient  was  the  same  for  both  bottom  layers,  the  energy 
ratios  were  summed  and  the  bottom-absorption  coefficients  were 
calculated  from  experimental  mode-attenuation  coefficients  as  described 
in  Chapter  II.  The  results  of  these  calculations  are  presented  in 
Table  11  for  the  two  short-range  legs  (G  and  H)  and  the  long-range 

leg  (J). 

4.4.4  Predicted  bottom  absorption — Site  2.  As  noted  in  the 
discussion  on  Site  1,  grab-sample  or  core  data  is  required  to  obtain 
input  parameters  for  the  sediment  model.  At  Site  2,  the  grab-samples 
noted  previously  (reported  by  Milliman  et  al.)  provided  adequate 
sampling,  particularly  near  the  long-range  propagation  leg.  In 


TABLE  10 


FREQUENCY  (Hz) 


Figure  20.  Daytona  Beach  Experimental  and 
Model  Dispersion — Site  2 
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TABLE  11 

ENERGY  DISTRIBUTION  (y)  AND  BOTTOM- 

ABSORPTION  COEFFICIENTS  (dj— SITE  2 

A 

Daytona  Beach 


Frequency 

(Hz) 

y(2> 

50 

0.07091 

0.02548 

63 

0.04297 

0.00592 

100 

0.02211 

0.00097 

160 

0.00909 

0.00008 

200 

0.00551 

0.00002 

250 

0.00319 

— 

315 

0.00175 

— 

400 

0.00091 

— 

500 

0.00048 

_  _ 

aA 

(dB/km) 

Y 

G 

H 

J 

0.09639 

26.97 

54.78 

— 

0.04889 

36.41 

101.66 

— 

0.02308 

55.46 

126.08 

— 

0.00917 

100.33 

121.05 

81.79 

0.00553 

144.67 

135.62 

117.54 

0.00319 

194.36 

244.51 

156.74 

0.00175 

285.71 

23*+. 29 

245.71 

0.00091 

439.56 

439.56 

417.58 

0.00048 

666.67 

666.67 

708.33 

addition,  a  single  core  obtained  by  the  Naval  Research  Laboratory  (NRL) 
is  indicated  near  Site  2  in  Figure  10.  By  selecting  those  grab-samples 
closest  to  the  point  where  Leg  J  crosses  the  sample  line,  the  data  in 
Table  12  was  obtained.  Data  for  the  NRL  core  (averaged  over  the  40  cm 
core  length)  is  also  presented  in  the  table.  From  the  average  mean 
grain-size  (Mz)  and  sorting  (0^)  for  the  grab-samples,  the  permeability 
(K)  and  pore-size  parameter  (a)  were  calculated  and  the  sediment  model 
input  parameters  in  Table  13  assumed.  The  bulk  modulus  for  the  mineral 
frame  (kr)  was  calculated  using  the  VRH  method  by  assuming  a  mineral 
content  of  60  percent  calcium  carbonate  and  40  percent  quartz.  With 
these  input  parameters,  the  Biot  sediment  model  was  run  and  the 
absorption  values  in  Table  14  were  predicted. 

4.4.5  Discussion  of  results — Site  2.  The  absorption  predicted 
by  the  sediment  models  and  that  obtained  from  experimental  mode- 
attenuation  coefficients  is  compared  in  Figure  21.  Also  shown  is  the 
absorption  predicted  by  Hamilton  for  the  sediment  mean  grain-size  found 
at  Site  2.  As  may  be  seen,  the  absorption  predicted  by  both  techniques 
is  somewhat  lower  than  that  obtained  experimentally  below  150  Hz.  As 
noted  in  the  discussion  on  Site  1,  this  is  probably  due  to  the  depth  of 
penetration  into  the  bottom  occurring  at  the  lower  frequencies  [note 
Y  in  Table  18].  The  volume  absorption  coefficients  for  the  deeper 
layers  are  unknown;  however,  it  is  probable  that  attenuation  due  to 
conversion  to  shear  waves  occurs  since  McDaniel  (1980)  found  that  the 
shear  velocity  of  the  second  bottom  layer  is  higher  (670  m/sec)  than 
the  value  assume  here.  At  frequencies  above  150  Hz,  the  energy  tends 
to  be  concentrated  in  the  first  bottom  layer  and,  as  may  be  seen  in  the 
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TABLE  12 

GRAB- SAMPLE  AND  CORE  DATA— SITE  2 
Daytona  Beach 


Sample 

Number 

M 

z 

(phi) 

°<}> 

(phi) 

Percent 

CaC03 

8 

P 

(g/cm3) 

1379 

0.43 

1.77 

62.1 

— 

-  - 

1380 

0.96 

0.65 

47.3 

— 

— 

1381 

1.20 

0.69 

56.2 

— 

— 

1382 

0.11 

1.76 

64.5 

— 

— 

1383 

1.33 

0.44 

62.3 

— 

— 

1384 

1.02 

0.59 

55.1 

— 

— 

1385 

1.C5 

0.47 

52.5 

-  - 

— 

1386 

1.15 

0.58 

59.3 

— 

-  - 

1387 

0.61 

0.58 

58.9 

— 

— 

1388 

0.66 

0.74 

66.4 

— 

— 

Average 

0.85 

0.83 

58.5 

— 

-  - 

NRL 

Core 

0.59 

1.21 

— 

0.38 

2.06 

Note : 

phi  =  -  log,d 

;  d  is  the 

mean  grain 

diameter  in  mm. 
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TABLE  13 

SEDIMENT  MODEL  PARAMETERS-- SITE  2 


Daytona  Beach 


z 

(phi) 

0.85 


v  2 

(phi)  (cm  ) 
0.83  7.5  x  10' 


6 

0.38 


p  r 

3  r  2 

(g/cm  )  (dyne /cm  ) 

2.06  5.6  x  IQ11 


Note:  Other  input  parameters  identical  to  Table  6. 


f 

TABLE  14 

PREDICTED  SEDIMENT  ABSORPTION  (aA)~SITE  2 
Daytona  Beach 

“A 

(dB/km) 

40  Hz  100  Hz  200  Hz  300  Hz  400  Hz  500  Hz 
0.85  0.83  6.6  39.3  141.8  278.2  421.1  553.5 


FREQUENCY  (Hz) 


Figure  21.  Comparison  of  Experimental  and 
Predicted  Bottom  Absorption — 
Daytona  Beach  Site  2 
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figure,  Che  predicted  and  experimental  absorption  agree  much  better. 

The  absorption  predicted  by  Hamilton's  regression  formula  matches 
somewhat  better  at  low  frequency  but  gives  poorer  agreement  than  the 
Biot  model  at  high  frequency. 

4.5  Comparison  of  Results  at  Site  1  and  Site  2 

Although  Sites  1  and  2  are  separated  by  only  about  37  km,  some 

difference  does  exist  in  the  surficial  sediments  which  appears  to  be 

reflected  in  the  modal  attenuation.  As  discussed  in  the  previous 

sections,  the  surficial  sediment  at  Site  1  is  best  characterized  as  a 

moderately-sorted  (0,  =  0.84),  medium  sand  (M  *  1.44)  with  a  calcium 

(p  z 

carbonate  content  of  approximately  30  percent.  A  coarser  sediment 
(M  -  0.86)  with  nearly  the  same  sorting,  but  with  a  calcium  carbonate 
content  twice  that  of  Site  1,  exists  at  Site  2.  Based  on  dispersion 
studies,  the  uppermost  layers  are  on  the  average  very  thin,  i.e.,  4  m 
at  Site  1  and  6  m  at  Site  2.  The  compressional  sound  speed  in  these 
layers  is  consistent  with  the  average  grain-sizes  in  that  the  finer 
sediment  has  the  lower  sound  speed.  In  both  cases,  the  next  layer  is 
relatively  thick;  however,  the  compressional  velocity  of  the  half-space 
assumed  for  Site  2  is  somewhat  larger  (2200  m/sec)  as  compared  to  that 
of  the  second  bottom  layer  at  Site  1  (1800  m/sec). 

The  mode-attenuation  coefficients  measured  at  Sites  1  and  2  are 
comparable  (except  for  Leg  H  at  Site  2)  at  low  frequency;  however,  at 
the  higher  frequencies,  the  mode  attenuation  at  Site  1  is  roughly  twice 
that  found  at  Site  2.  Because  of  this.  Leg  H  at  Site  2  and  the  750  Hz 
data  at  Site  1  have  been  deleted  from  further  comparisons. 
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Kibblewhite  and  Denham  (1966)  note  that,  at  high  frequencies, 
the  mode  attenuation  for  the  first  mode  6^  should  follow  the  relation 


x  „  “A 

■L  H 


(49) 


Thus,  for  the  same  sediment  type,  the  mode  attenuation  Is  proportional 

_3 

to  H  ,  where  H  is  the  water  depth.  It  follows  then  that  the  mode 
attenuation  at  Site  1  should  be  larger.  If  the  ratio  of  the  mode 
attenuation  at  Site  1  to  that  at  Site  2  is  calculated  based  on 
Equation  (49)  (assuming  or^  the  same) ,  then  the  attenuation  at  Site  1 
should  be  3.8  times  that  at  Site  2.  An  examination  of  Table  3  and 
Table  9  indicates  that  this  does  not  hold  true;  hence,  if  we  assume  that 
Equation  (49)  is  valid,  then  the  bottom-absorption  coefficients  (aA)  for 
the  two  sites  must  be  different.  This  is  confirmed  by  the  bottom- 
absorption  coefficients  calculated  from  the  mode-attenuation  coeffici¬ 
ents.  These  differences  are  readily  apparent  in  Table  15  where  power- 
law  fits  have  been  calculated  for  the  experimental  and  the  predicted 
bottom-absorption  coefficients.  The  experimental  relations  for  Site  2 
show  some  variation  in  magnitude.  This  may  be  due  to  the  high  shear 
speed  recently  determined  for  this  leg  by  McDaniel.  The  differences 
noted  in  sediment  type  for  the  two  sites,  although  small,  are  reflected 
in  the  bottom  absorption  predicted  by  both  the  Biot  and  Hamilton 
methods.  Over  the  frequency  ranges  discussed  (50  to  600  Hz  at  Site  1 
and  40  to  500  Hz  at  Site  2) ,  the  absorption  predicted  by  the  Biot 
theory  provides  better  agreement  over  the  greatest  proportion  of  the 
frequency  range.  This  appears  to  indicate  that  a  frequency  dependence 
other  than  f^  is  required. 


I 
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Site 

1 

2 


TABLE  15 


POWER-LAW  FITS  TO  BOTTOM-ABSORPTION 

COEFFICIENTS  (a.)— SITES  1  AND  2 
A 

Daytona  Beach 


Experimental  a 

iV 

Leg  (dB/km) 

E  0.013  f1'65 

F  0.013  f1,65 


Predicted 

(dB/km) 

Biot 

0.004  f1,84 


G  0.179  f1'28  0.011  f1'76 

J  0.005  f1*88 


Hamilton 
0.49  f 

0.47  f 


Note:  The  frequency  f  is  in  Hertz 


CHAPTER  V 


BOTTOM- ABSORPTION  ANALYSIS 
CORPUS  CHRISTI  EXPERIMENTS 


5.1  Introduction 

In  the  previous  chapter,  data  were  analyzed  from  a  test  site 
where  the  sediments  were  In  the  medium-to-coarse  sand  range.  The  clay 
and  slit  sediments  found  off  the  coast  of  Corpus  Christi,  Texas,  provide 
test  data  at  the  extreme  end  of  the  grain-size  scale.  Cores  were 
obtained  at  the  test  site  and  were  analyzed  for  sediment  parameters 
which  provide  the  basis  for  bottom-absorption  predictions. 

The  procedure  followed  is  much  the  same  as  that  followed  in  the 
preceding  chapter.  The  experimental  program  is  defined,  and  received 
energy  and  mode-attenuation  coefficients  calculated.  The  bottom  model 
derived  is  somewhat  unique  for  continental-shelf  areas  in  that  the 
uppermost  sedimentary  layer  has  a  sound  speed  less  than  that  of  the 
water  column.  This  model  is  used  in  the  normal-mode  program  to  calculate 
energy  ratios  and,  hence,  experimental  bottom  absorption.  The  experi¬ 
mental  values  are  then  compared  with  absorption  predicted  by  the 
sediment  models. 

5.2  Experiments 

During  April  1977,  experiments  were  performed  off  the  coast  of 
Corpus  Christi,  Texas.  Short-range  and  long-range  propagation  and 
seismic-refraction  experiments  were  performed  as  described  in  Chapter 
III;  however,  only  constant-depth  leg  results  are  considered  here. 


Both  short-range  and  long-range  propagation  data  were  analyzed 
by  Rubano  (1980)  for  the  site  and  legs  shown  In  Figure  22.  The  water 
depth  at  the  receiver  was  30.5  m  and  varied  by  less  than  0.5  m  along 
the  short-range  constant-depth  leg.  The  maximum  range  attained  was 
14.7  km.  For  the  short-range  tests,  small  (7.1  g)  charges  were 
detonated  at  depths  of  9.5  and  19  m.  The  new  vertical  array  was  used 
as  the  receiver. 

Long-range  data  were  obtained  using  an  aircraft  to  drop  SUS 
charges,  set  to  detonate  at  a  depth  of  18.3  m,  out  to  a  range  of  77.6  km 
along  the  constant-depth  track.  The  bottom-sensor  unit  was  in  place  for 
the  long-range  tests;  hence,  the  top  hydrophone  (1  m  above  the  bottom) 
was  used  as  the  receiver  for  these  tests.  Along  this  leg,  the  bathymetry 
varied  as  shown  in  Figure  23. 

Seismic-refraction  analysis  was  performed  on  data  obtained  along 
the  constant-depth  leg  out  to  12  km  in  range.  A  12  kHz  fathometer  was 
towed  over  the  same  leg,  yielding  bottom  and  shallow  sub-bottom  data. 

A  total  of  four  cores,  which  will  be  discussed  later,  were  obtained  for 
this  site.  Sound-velocity  profiles  obtained  along  both  tracks  showed  a 
very  slight  positive  gradient. 


5.3  Analysis  of  Data 

The  received  energy  (corrected  for  cylindrical  spreading)  was 
calculated  for  short-range  and  long-range  data,  and  mode-attenuation 


coefficients  were  obtained  for  the  first  normal  mode.  The  bottom  model 


derived  by  Rubano  (1980),  from  dispersion  and  seismic-refraction 
analysis,  was  used  in  a  four-layer  normal-mode  program  to  calculate 
the  energy  ratios  for  each  layer.  Experimental  bottom-absorption 


Figure  22.  Corpus  Christ!  Test  Area 


DEPTH  (ml 
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coefficients  were  then  calculated.  Bottom  absorption  was  predicted  by 
the  sediment  models,  using  input  data  based  on  the  analysis  of  the  cores 
obtained  at  the  site,  and  compared  with  the  experimental  results. 

5.3.1  Received-energy  and  mode-attenuation  calculations. 
Short-range  data  were  obtained  for  ranges  from  1.0  to  14.7  km  and  for 
frequencies  from  25  to  400  Hz.  Energy  received  on  the  bottom  hydrophone 
of  the  array  was  analyzed  for  the  small  charges  detonated  at  a  depth  of 
19.5  m.  Rubano  (1980)  noted  that  the  portion  of  the  signal  due  to  the 
first  mode  was  well  separated  in  time  from  the  higher-order  mode 
arrivals  (basically,  the  second  mode).  Hence,  the  data  from  a  single 
hydrophone  of  the  array  was  processed  using  temporal  separation  to 
obtain  the  first  mode.  Two  shots  were  analyzed  for  each  of  ten  range 
points  (at  most  frequencies)  spaced  over  the  range  interval  noted  above. 
Both  shots  were  averaged  to  obtain  the  plot  of  relative  received  energy 
shown  in  Figure  24.  A  least-squares,  straight-line  fit  was  calculated 
for  all  of  the  data  points  at  each  frequency  and  is  also  plotted  in  the 
figure.  The  mode-attenuation  coefficients  obtained  from  the  slope  of 
this  straight  line  are  listed  in  Table  16  with  the  calculated  95  percent 
confidence-interval  limits  for  the  slope.  Note  that  only  data  up  to 
250  Hz  is  presented,  whereas  received  energy  was  calculated  out  to 
400  Hz.  This  limitation  is  imposed  due  to  unusual  propagation  effects, 
which  will  be  discussed  later,  that  occurred  at  the  higher  frequencies. 

Long-range  data  were  analyzed  for  ranges  from  18.5  to  55.6  km 
over  a  nominal  frequency  range  of  25  to  250  Hz.  Shots  were  detonated 
to  a  range  of  77.6  km;  however,  beyond  55.6  km,  the  signal-to-noise 
ratio  was  inadequate  for  data  analysis.  The  relative  received  energy 


104 


TABLE  16 

SHORT-RANGE  MODE-ATTENUATION  COEFFICIENT  (6) 
Corpus  Christi 


Mode-Attenuation 

95%  Confidence 

Interval 

Frequency 

(Hz) 

Coefficient 

(dB/km) 

Upper  Limit 

Lower  Limit 

25 

0.48 

0.64 

0.32 

31.5 

0.45 

0.58 

0.32 

40 

0.38 

0.50 

0.26 

50 

0.36 

0.49 

0.24 

63 

0.33 

0.45 

0.22 

80 

0.45 

0.59 

0.32 

100 

0.56 

0.70 

0.42 

125 

0.53 

0.71 

0.35 

160 

0.83 

1.03 

0.63 

200 

1.13 

1.30 

0.96 

250 


1.29 


1.42 


1.15 
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was  calculated  using  corrections  for  gain  and  cylindrical  spreading 
and  is  plotted  in  Figure  25  for  each  frequency.  The  mode-attenuation 
coefficients  obtained  from  the  straight-line  fits  to  these  data  are 
listed  in  Table  17. 

5.3.2  Bottom  model.  The  bottom  model  derived  by  Rubano  is 
presented  in  Figure  26.  The  four-layer  model  is  based  on  seismic- 
refraction  and  dispersion  analysis.  Note  that  the  upper  sediment  layer 
has  a  sound  speed  less  than  that  found  in  the  water  column.  This  low 
velocity  (1490  m/sec)  is  consistent  with  the  grain-size  and  velocity 
measurements  obtained  from  analysis  of  the  four  cores.  The  effect  of 
the  low  sound  speed  in  this  layer  is  to  trap  bottom-incident  energy 
above  a  frequency  (estimated  by  Rubano)  of  250  Hz.  This  is  recognized 
as  the  angle  of  intromission  effect  associated  with  the  plane-wave 
reflection  and  transmission  coefficients,  in  which  the  reflected  signal 
goes  to  zero  for  certain  frequencies  and  angles  of  incidence.  Since  the 
measured  losses  are  directly  related  to  this  phenomenon,  the  frequency 
range  for  the  data  presented  is  limited  to  250  Hz. 

The  bottom  model  derived  by  Rubano  appears  to  be  consistent  with 
other  data  from  the  literature.  Although  the  closest  seismic-refraction 
data  to  be  found  are  roughly  300  km  to  the  northeast,  near  Galveston, 
Texas,  the  velocities  shown  in  Table  18  are  in  general  agreement  with 
the  model.  These  studies  by  Ewing,  Antoine,  and  Ewing  (1960),  and 
Antoine  and  Ewing(1963),  also  agree  in  general  with  seismic-refraction 
results  obtained  very  near  the  test  site  in  water  18  m  deep.  A  low- 
velocity  layer  is  not  reported  in  these  references.  However,  seismic- 
refraction  analysis  assumes  that  the  velocity  always  increases  with 
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TABLE  17 


LONG-RANGE  MODE-ATTENUATION  COEFFICIENT  (5) 
Corpus  Christ! 


Frequency 

(Hz) 

Mode-Attenuation 

Coefficient 

(dB/km) 

95%  Confidence 

Upper  Limit 

Interval 

Lower  Limit 

25 

0.86 

1.09 

0.63 

31.5 

0.56 

0.65 

0.47 

40 

0.47 

0.51 

0.42 

50 

0.43 

0.49 

0.38 

63 

0.40 

0.43 

0.36 

80 

0.46 

0.51 

0.42 

100 

0.56 

0.61 

0.51 

125 

0.67 

0.73 

0.62 

160 

0.72 

0.79 

0.65 

200 

0.73 

0.83 

0.63 

250 


0.79 


0.86 


0.72 
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depth;  hence,  the  referenced  refraction  results  do  not  include  the 
low-velocity  layer.  In  addition,  layer  thicknesses  will  be  in  error 
due  to  neglect  of  the  low-velocity  layer.  Thus,  the  bottom  model  of 
Rubano  is  considered  to  be  valid. 

5.3.3  Experimental  bottom  absorption.  Using  the  four-layer 
bottom  model  and  normal-mode  program,  energy  ratios  for  the  sub-bottom 
layers  were  calculated.  These  values  are  shown  in  Table  19  where  Y  is 
the  sum  of  the  values  for  the  individual  layers.  With  <5  known  for 
each  frequency,  the  experimental  bottom  absorption  (a  )  was  calculated 
by  dividing  the  mode-attenuation  coefficients  for  the  first  mode  in 
Tables  16  and  17  by  the  appropriate  value  of  Y  .  These  coefficients 
are  also  listed  in  Table  19. 

5.3.4  Predicted  bottom  absorption.  During  the  experimental 
program,  four  cores  were  obtained  near  the  30  m  test  site.  The  cores 
were  analyzed  by  the  Geological  Laboratory  of  the  Naval  Oceanographic 
Office  where  velocities  ranging  from  1490  to  1500  m/sec  were  measured 
through  the  cores.  Grain-size  measurements  on  the  cores  show  that  this 
thin  layer  consists  of  a  very  poorly  sorted,  medium-to-f ine  silt.  A 
summary  of  the  data  provided  by  the  core  analysis  is  listed  in  Table  20. 
The  grain-sizes  listed  represent  the  maximum  and  minimum  values  within 
each  core  except  for  core  number  10  where  only  one  value  was  used.  The 
dominant  constituents  of  the  cores  were  quartz  sand  and  shell  fragments. 
Cores  obtained  by  other  researchers  are  indicated  in  Figure  22.  Trask 
(1953)  analyzed  these  cores  for  organic  content,  chemical  composition, 
and  water  content,  and  found  median  grain-sizes  averaging  7.5  phi  and 
calcium  carbonate  content  averaging  3.4  percent  near  the  numbered  cores. 
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TABLE  19 

ENERGY  DISTRIBUTION  (y)  AND 

BOTTOM-ABSORPTION  COEFFICIENTS  (a.) 

A 

Corpus  Christ! 


(dB/km) 


Frequency 

(Hz) 

Y(1> 

(2) 

y 

(3) 

Y 

y 

Short 

Range 

0.83 

Long 

Range 

1.48 

25 

0.1016 

0.3886 

0.0897 

0.5799 

31.5 

0.1145 

0.2816 

0.2230 

0.4184 

1.08 

1.34 

AO 

0.1143 

0.1788 

0.0043 

0.2974 

1.28 

1.58 

50 

0.1074 

0.1119 

0.0007 

0.2200 

1.64 

1.95 

63 

0.0974 

0.0681 

8.2  x  10-5 

0.1656 

1.99 

2.42 

80 

0.0868 

0.0410 

5.9  x  io"6 

0.1278 

3.52 

3.60 

100 

0.0780 

0.0258 

— 

0.1038 

5.39 

5.39 

125 

0.0711 

0.0166 

— 

0.0877 

6.04 

7.64 

160 

0.0668 

0.0106 

— 

0.0774 

10.72 

9.30 

200 

0.0673 

0.0076 

— 

0.0749 

15.09 

9.75 

250 

0.0756 

0.0062 

—  - 

0.0818 

15.77 

9.66 
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Core 

Number 


TABLE  20 

CORE  DATA-MAXIKUM  AND  MINIMUM  GRAIN  SIZES 
Corpus  Christ! 

0, 


Length 

(cm) 


M 

z 

(phi) 


$ 

(phi) 


P‘3 
(g/cin  ) 


P  3  B 

(g/cin  )  _ 


6 

47 

6.07 

2.00 

2.67 

1.82 

0.51 

6 

47 

8.78 

2.13 

2.73 

1.52 

0.70 

9 

69 

6.37 

2.19 

2.67 

1.58 

0.65 

9 

69 

8.78 

2.10 

2.73 

1.48 

0.72 

10 

59 

8.67 

2.17 

2.71 

1.51 

0.70 

15 

36 

6.51 

2.61 

2.67 

1.73 

0.56 

15 

36 

8.38 

2.36 

2.68 

1.57 

0.66 

These  results  suggest  that  the  data  In  Table  20  is  representative  of 
the  test  area  and  that  the  mineral  composition  of  the  sediment  frame 
may  be  considered  pure  quartz  (since  the  carbonate  content  is  low) . 

Using  the  data  in  Table  21,  sediment  model  runs  (with  the  Biot 
model)  were  made  to  predict  the  absorption  within  the  sedimentary  layer. 
The  permeability  (K)  and  pore-size  parameters  (a)  were  calculated  from 
the  mean  grain-size  (M  )  and  the  standard  deviation  (a.)  as  described 

Z  m 

in  Chapter  II.  The  densities  and  porosities  are  measured  data  from  the 
core  analysis.  Note  that  values  different  than  those  used  in  Chapter  IV 
are  used  for  a  ,  Ag  and  Ag  for  the  fine  sediments  found  at  this 
site.  These  values  are  recommended  by  Stoll  (1974)  as  providing  the 
best  fit  to  the  fine-grained  experimental  data  of  Wood  and  Weston  (1964). 

The  results  of  the  sediment  model  runs  are  listed  in  Table  21. 

As  may  be  seen,  the  range  of  absorption  values  predicted  is  rather 
small  for  the  mean  grain-sizes  used.  If  the  absorption  is  expressed 
in  the  power  law  form,  fn  ,  by  calculating  a  power-law  fit  to 

the  absorption  as  a  function  of  frequency,  the  exponent  n  is  one  in 
all  cases  and  KA  varies  between  0.052  and  0.061  dB/km*Hz.  A  typical 
value  for  the  absorption  is  aA  =  0.057  f  .  Rubano  found  by  trial  and 
error  that  the  bottom  absorption  was  0.005  f ^ ^  nepers/m  or  in  the 
units  used  here,  aA  *  0.022  f^*^  dB/km  . 

5.4  Discussion  of  Results 

The  results  of  the  experimental  and  predicted  bottom-absorption 
calculations  are  plotted  in  Figure  27.  Two  different  prediction 
equations,  from  the  models  of  Hamilton  and  Biot,  are  compared  with  both 
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Figure  27.  Comparison  of  Experimental  and  Predicted 
Bottom  Absorption — Corpus  Christi 
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short-range  and  long-range  calculations  of  bottom  absorption  from 
Table  19.  Both  model  predictions,  shown  by  the  solid  lines,  are  based 
on  the  largest  mean  grain-size  for  core  9  (6.37  phi).  The  range  of 
values  predicted  by  the  Biot  model,  based  on  the  mean  grain-sizes  from 
the  core  data* is  indicated  by  the  dotted  lines  in  the  figure.  If  a 
prediction  based  upon  porosity  is  calculated  with  the  Hamilton  model, 
the  values  are  even  larger.  As  can  be  seen,  the  predictions  by  Hamilton 
are  much  too  large  over  the  entire  frequency  range,  whereas  the  predic¬ 
tions  based  on  the  Biot  model  give  good  agreement  except  at  low 
trequencies.  It  might  be  expected  that  the  predictions  are  poorest 
at  the  low  frequencies  where  significant  energy  penetrates  to  the 
second  and  third  layers.  In  calculating  the  experimental  bottom 
absorption,  it  was  assumed  that  the  attenuation  was  the  same  in  all  of 
the  bottom  layers.  This  is  probably  not  true  and,  at  low  frequencies 
with  deep  penetration,  this  is  reflected  in  poorer  agreement  between 
experimental  and  predicted  absorption.  For  the  data  in  Figure  27,  the 
best  agreement  occurs  from  100  to  250  Hz  for  the  short-range  data. 

Over  this  frequency  range,  75  to  92  percent  of  the  bottom  energy  is  in 
the  first  layer,  as  is  shown  in  Table  19.  Thus,  it  appears  that  the 
predictive  ability  of  the  sediment  models  is  best  assessed  over  this 
frequency  range  at  the  Corpus  Christi  test  site. 


CHAPTER  VI 


BOTTOM-ABSORPTION  ANALYSIS 
PANAMA  CITY  EXPERIMENTS 


6.1  Introduction 

Using  mode-attenuation  data  from  the  literature,  experimental 
and  predicted  bottom  absorption  was  calculated  for  constant-depth 
propagation  paths  in  water  depths  of  18.3  and  30.5  m.  The  seabed  in 
Panama  City  test  areas  consisted  of  moderately  well-sorted,  fine-to- 
medium  sand  with  a  carbonate  content  of  approximately  18  percent.  The 
area  was  similar  to  the  Daytona  Beach  test  area  in  terms  of  the 
surficial  sediments,  but  differed  in  the  sub-bottom  structure.  A 
separate  environmental  measurement  and  analysis  program  conducted  by 
the  University  of  Texas  at  Galveston  provided  new  and  refined  data 
which  were  used  in  modeling  the  sub-bottom  structure. 

6.2  Experiments 

Shallow-water  propagation  experiments  were  performed  off  the 
coast  of  Panama  City,  Florida,  in  January,  April,  and  June  of  1976. 

The  experiments  to  be  discussed  here  were  performed  using  the  vertical 
array  (LC-10  array)  during  the  June  time-period.  Mode-attenuation 
data  were  obtained  along  the  two  short-range  propagation  legs  shown 
in  Figure  28.  Seismic-refraction  data  were  also  obtained  and  analyzed. 
Cores  were  obtained  and  analyzed  by  the  University  of  Texas  at 


Galveston. 


Figure  28.  Panama  City  Test 
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At  each  site,  the  receiving  array  was  cabled  directly  to  an 
oceanographic  platform:  Stage  I  in  30.5  m  of  water,  and  Stage  II  in 

18.3  m  of  water.  From  Stage  I,  short-range  data  were  obtained  out  to 
a  maximum  range  of  17.8  km  along  a  constant-depth  track  (Leg  C)  in 
water  30.5  m  deep.  Short-range  data  were  obtained  from  Stage  II  out 
to  a  maximum  range  of  17  km  along  constant-depth  Leg  D  (18.3  m  water 
depth).  The  receiving  array  for  the  experiment  from  Stage  II  was 
shortened,  using  five  of  the  eight  hydrophones;  the  remaining  three 
were  coiled  at  the  surface.  The  signal  sources  for  both  experiments 
were  the  small  (7.1  g)  charges.  These  charges  were  detonated  at  depths 
of  9.5  m  and  19.2  m  along  Leg  C,  and  5.8  and  11.3  m  along  Leg  D.  In 
both  cases,  the  shallow  source  was  predicted  to  (and,  in  fact,  did) 
maximize  the  excitation  of  the  second  mode,  while  the  deeper  source 
was  designed  to  minimize  the  second  mode. 

The  water  depth  along  both  tracks  was  determined  to  be  essential¬ 
ly  constant.  The  sound-velocity  profiles  obtained  during  the  experiment 
showed  only  a  slight  variation  with  range.  The  profiles  used  in  the 
analysis  were  typical  of  those  obtained  and  will  be  presented  with  the 
bottom  models  in  a  later  section. 

6.3  Analysis  of  Data 

Propagation  along  both  Legs  C  and  D  was  analyzed  extensively  by 
Caswell  (1979).  Using  array  processing  of  the  data  to  select  the 
desired  modes,  he  obtained  mode-attenuation  coefficients  for  the  first 
mode  over  a  frequency  range  of  50  to  750  Hz  along  Leg  C  and  100  to 
750  Hz  along  Leg  D.  The  attenuation  of  the  second  and  third  modes  was 
also  obtained  for  portions  of  the  data. 
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6.3.1  Mode-attenuation  coefficients.  By  processing  the 
short-range  data  along  Legs  C  and  D  using  the  array  techniques 
described  in  Chapter  III,  the  mode-attenuation  coefficients  listed 

in  Table  22  were  obtained  by  Caswell.  Along  Leg  C,  these  results  may 
be  compared  with  those  of  Ingenito  (1973).  At  400  Hz,  Ingenito  reported 
a  mode-attenuation  coefficient  of  0.61  +  0.19  dB/km  within  95  percent 
confidence  limits,  giving  reasonable  agreement  with  the  results  of 
Caswell  (0.85  +  0.08).  However,  at  750  Hz,  Ingenito  found  mode 
attenuation  of  0.30  +  0.39  dB/km,  while  Caswell  found  attenuation  of 
1.28  +  0.24  dB/km.  The  reason  for  the  difference  at  750  Hz  is  not 
apparent  from  either  paper. 

6.3.2  Bottom  model.  A  valid  bottom  model  is  necessary  to 
calculate  the  experimental  bottom-absorption  coefficients  from  the 
mode-attenuation  values  listed  in  Table  22.  An  incorrect  choir'-  jf 
bottom  parameters  can  lead  to  large  errors  in  the  calculate  .osorption 
coefficients  since  the  modal  energy  distribution  and,  hence,  the  energy 
ratios  are  dependent  upon  these  parameters.  The  bottom  models  selected 
by  Caswell  and  Ingenito  are  quite  different  as  may  be  seen  in  Table  23. 
The  key  parameter  here  is  the  sound  speed  of  the  first  bottom  layer 
(c^) .  While  substantial  differences  exist  in  the  choice  of  these 
parameters,  the  effect  of  these  differences  is  not  easily  assessed 
because  of  the  differences  in  mode  attenuation  noted  in  the  previous 
paragraph.  As- will  be  discussed  in  the  summary  of  this  chapter,  the 
results  of  bottom-absorption  calculations  by  Caswell  and  Ingenito 

are  remarkably  close  in  their  frequency  dependence. 
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TABLE  22 

MODE-ATTENUATION  COEFFICIENT  (6) 

[Caswell  (1979)] 

Panama  City 

Mode-Attenuation  95%  Confidence  Interval 


Frequency 

(Hz) 

Coefficient 

(dB/km) 

Upper  Limit 

Lower  Limit 

50 

1.80 

2.14 

1.46 

100 

1.22 

1.36 

1.08 

200 

1.01 

1.23 

0.89 

400 

0.85 

0.93 

0.77 

600 

0.91 

1.05 

0.77 

750 

1.28 

1.04 

1.52 

100 

2.25 

2.56 

1.94 

200 

1.96 

2.18 

1.74 

400 

1.70 

1.86 

1.54 

600 

1.77 

1.55 

1.99 

750 


1.68 


1.96 


1.40 
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TABLE  23 

OTHER  BOTTOM  MODELS  FOR  PANAMA  CITY 


Bottom  Parameters 


Leg 


1 

(m/sec) 


H1 

(m) 


Bottom 

Model 

Source 


C 


1711  3.8  2100 


Caswell  (1979) 


D  1708  ® 


Caswell  (1979) 


C  1589 


Ingenito  (1973) 
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The  bottom  parameters  in  Table  23  were  based  on  the  best 
available  information  at  the  time  the  analysis  was  performed.  More 
recently,  however,  results  of  refraction  and  reflection  analysis,  core 
analysis,  and  in-situ  sound-velocity  measurements  performed  by  personnel 
of  the  Marine  Science  Laboratory  of  the  University  of  Texas  at  Galveston 
have  become  available.  Addy  et  al.  (1979)  correlated  the  results  of 
sub-bottom  reflection  surveys  with  sediment  lithology  based  on  long 
piston  cores.  McMillan  et  al.  (1978)  discussed  the  sediments  of  the 
test  area  (again,  based  on  the  cores)  in  detail.  Velocity  measurements 
were  made  in-situ  using  a  velocimeter  placed  in  the  cutting  head  of  the 
coring  apparatus  [see  Shirley  and  Anderson  (1975)].  Using  data  from  the 
reflection  surveys,  Houston  (1978)  performed  a  refraction  analysis  of 
this  data  and  obtained  both  velocities  and  layering  of  the  sub-bottom 
structure  within  the  test  area. 

The  environmental  models  used  in  this  study  for  Legs  C  and  D  are 
listed  in  Table  24.  The  sound-velocity  profiles  are  those  used  by 
Caswell  as  representative  of  the  complete  propagation  path.  The  bottom 
models  differ  from  those  of  both  Caswell  and  Ingenito,  being  based  on 
the  work  noted  in  the  previous  paragraph.  The  first  layer  velocity  for 
Leg  C  (1640  m/sec)  is  substantiated  by  an  in-situ  velocity  measurement 
made  near  Stage  1.  The  thickness  of  this  layer  is  that  used  by  Caswell 
and  is  based  on  seismic-refraction  analysis.  This  layer  thickness 
appears  reasonable  based  on  the  length  of  the  cores  which  yielded  the 
velocity  (2  to  4  m) .  In  addition,  Houston  noted  a  layer  thickness  of 
15  m  for  a  layer  of  comparable  velocity  at  a  distance  of  33  km  along  a 
252°  bearing  from  the  stage.  The  second  velocity  (1850  m/sec)  is  an 
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TABLE  24 

ENVIRONMENTAL  MODELS — LEGS  C  AND  D 
Panama  City 


Leg  C 

Leg  D 

Depth 

(m) 

Sound 

Velocity 

(m/sec) 

Density 

Depth 

(m) 

Sound 

Velocity 

(m/sec) 

Density 

0 

1527.5 

0 

1527.0 

1.1 

1526.0 

1.0 

4.1 

1526.3 

1.0 

30.5 

1525.0 

10.5 

1526.2 

18.3 

1525.2 

- 

-  -  Water/Bottom  Interface  -  - 

- 

30.5 

1640.0 

18.3 

1640.0 

34.3 

1640.0 

2.0 

22.1 

1640.0 

2.0 

34.3 

1850.0 

22.1 

1850.0 

oo 

1850.0 

OO 

1850.0 
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average  of  those  along  the  reflection  leg  (Figure  28)  as  determined 
by  Houston. 

The  recent  studies  noted  above  covered  an  area  seaward  from 
Stage  I.  Hence,  these  findings  must  be  extrapolated  inshore  to 
establish  a  bottom  model  for  the  region  near  Stage  II  and  along  Leg  D. 
Ludwick  (1964)  studied  the  sediments  in  the  northeastern  Gulf  of  Mexico 
and  classified  the  sediment,  for  a  large  region  which  included  Panama 
City,  as  quartz  sand  of  thickness  less  than  15  m.  He  noted  the  seaward 
distribution  of  this  sand  alternated  between  fine  and  medium  sand  out 
to  a  depth  of  about  20  fathoms.  Calculations  based  on  the  grain-size 
distribution  of  Ludwick  for  sands  in  this  region  yielded  a  mean  grain- 
size  of  2.5  phi  and  a  standard  deviation  of  0.51  phi.  Gould  and 
Stewart  (1955)  also  studied  the  sediments  off  Northern  Florida  and 
arrived  at  basically  the  same  conclusion;  that  is,  that  the  sediments 
are  quartz  sands  with  no  progressive  change  in  particle-size  distribu¬ 
tion  with  increasing  distance  from  shore.  On  this  basis,  it  was 
concluded  that  the  bottom  model  for  Leg  D  should  not  differ  substanti¬ 
ally  from  that  of  Leg  C.  Thus,  the  same  basic  bottom  model  (see 
Table  24)  was  used  for  both  sites. 

6.3.3  Experimental  bottom  absorption.  With  bottom  models 
established  for  Legs  C  and  D,  it  is  now  possible  to  calculate  energy 
ratios  (7)  using  the  three-layer  normal-mode  model,  and  hence,  the 
experimental  bottom  absorption  (aA) .  These  calculated  values  are 


listed  in  Table  25  for  both  legs  using  the  mode-attenuation  coefficients 
from  Table  22  measured  by  Caswell.  If  the  mode-attenuation  coefficients 
of  Ingenito  are  used  and  the  bottom  absorption  is  calculated  with 


TABLE  25 
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the  energy  ratios  from  Table  25,  bottom-absorption  values  of 
260.7  dB/km  at  400  Hz  and  489.4  dB/km  at  750  Hz  are  obtained. 

6.3.4  Predicted  bottom  absorption.  As  noted  previously,  the 
sediments  in  the  vicinity  of  Stage  I  have  been  described  in  some  detail 
by  Addy  et  al. (1979),  and  McMillan  et  al.  (1978).  Both  studies  are 
based  on  piston  cores  obtained  along  sample  lines  emanating  from  Stage  I. 
Within  8  km  of  Stage  I,  four  cores  were  obtained  with  lengths  ranging 
from  2  to  4  m.  The  sediments  within  the  cores  were  generally  classified 
by  McMillan  as  quartz-shell  sands  having  mean  grain-sizes  ranging  from 
1.35  to  1.85  phi,  with  an  average  of  1.64.  The  sorting  of  this  type  of 
sand  ranged  from  a  standard  deviation  of  0.53  to  0.80  phi  and  the 
calcium  carbonate  content  averaged  approximately  18  percent.  Addy 
estimated  the  mean  grain-size  for  a  larger  area  as  2.3  phi  with  a 
standard  deviation  of  1.0  phi.  Using  the  grain-size  values  of  McMillan, 
the  Biot  sediment  model  was  used  to  predict  the  bottom  absorption  for 
frequencies  from  50  to  750  Hz.  The  grain-sizes,  other  bottom  parameters, 
and  predicted  absorption  values  are  listed  in  Table  26  (a)  and  (b) . 

Since  some  evidence  exists  for  finer  grain  sizes  in  the  test  area,  a 
mean  grain-size  of  1.85  phi  [from  Table  26  (a)]  was  used  as  being 
typical  for  the  area. 

6.4  Discussion  of  Results 

Mode-attenuation  data  from  papers  by  Caswell  and  Ingenito  were 
used  for  this  test  area  to  obtain  experimental  bottom-absorption  data. 
Recent  studies  of  the  sediments  off  Panama  City  indicated  that  the 
bottom  models  assumed  by  both  researchers  might  be  improved.  Based  on 
these  studies,  a  new  bottom  model  was  derived  for  both  Legs  C  and  D, 


TABLE  26  (a) 


BIOT  SEDIMENT  MODEL  PARAMETERS 
Panama  City 


Bottom  Parameters 


Sample 

No. 

M 

z 

(Phi) 

°4> 

(phi) 

K 

(cm2) 

a 

(cm) 

8 

P 

(g/cm 

1 

1.35 

0.80 

4.1  x  lo"7 

2.5  x  10'3 

0.375 

2.07 

2 

1.64 

0.62 

3.5  x  10“7 

2.3  x  10“3 

0.400 

2.03 

3 

1.85 

0.53 

2.9  x  io-7 

2.1  x  IQ-3 

0.410 

2.01 

11  2 

Note:  =  4.2  x  10  dyne/cm 

10  2 

*  2.4  x  10  dyne/cm 

Pr  =*  2.71  g/cra3 


-2  2 
H  =  1.0  x  10  dyne-sec/cm 


v  ®  3  x  10  cm/sec  A„  =  0.15 

E  E 

v  *  2.1  *  104  cm/sec  A  *  0.20 
s  s 


1.00  g/ cm' 


a  *  1.25 


TABLE  26  (b) 

PREDICTED  BOTTOM  ABSORPTION  (a.) 

A 

Panama  City 


_ Absorption  (dB/km) 

Sample 


No. 

50  Hz 

100  Hz 

200  Hz 

400  Hz 

600  Hz 

750  Hz 

1 

6.1 

22.4 

83.5 

294.6 

563.6 

768.1 

2 

5.2 

18.6 

69.5 

252.6 

507.2 

717.8 

3 

4.4 

15.7 

58.2 

214.8 

444.0 

644.4 
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and  the  experimental  bottom  absorption  was  then  calculated  for 
frequencies  from  50  to  750  Hz.  The  most  sensitive  parameter  in  the 
bottom  model  was  determined  to  be  the  velocity  of  the  first  bottom 
layer.  For  example,  as  this  parameter  was  varied  from  1620  to  1711 
m/sec,  the  absorption  ( a varied  from  0.022  f^"83  to  0.007  f3"^  dB/km 
(f  in  Hz).  On  the  other  hand,  as  the  first  layer  thickness  varied  from 
3.8  to  10  m,  the  absorption  CctA>  varied  from  0.015  f1'73  to  0.0065  f1'86 
dB/km.  The  bottom  model  chosen  for  Legs  C  and  D  (Table  24)  used  a 
sound  speed  of  1640  m/sec  determined  from  in-situ  velocity  measurements 
on  cores  near  Stage  I.  Using  this  bottom  model  and  the  mode-attenuation 
coefficients  of  Caswell  and  Ingenito,  the  bottom  absorption  was 
calculated.  From  analysis  of  these  same  cores,  mean  grain-size  data 
were  obtained  and  used  in  the  sediment  models  of  Biot  and  Hamilton  to 
predict  bottom  absorption. 

The  experimental  and  predicted  bottom-absorption  data  is  plotted 
in  Figure  29  and  listed  in  Table  27  in  the  power  law  form  a  *  K  fn  . 

A  A 

Also  listed  in  Table  27  is  the  original  data  of  Caswell  and  Ingenito. 

In  the  figure,  the  Hamilton  model  matches  the  experimental  absorption 
reasonably  well  at  low  frequencies,  whereas  the  Biot  model  does  not 
predict  enough  absorption.  At  the  higher  frequencies  (600  to  750  Hz), 
neither  model  predicts  enough  absorption  when  compared  with  the  results 
obtained  from  Caswell's  mode-attenuation  data.  When  the  predicted 
values  are  compared  with  the  experimental  values  based  on  Ingenito 's 
measurements,  the  agreement  is  better.  Because  of  this  inconsistency, 
the  750  Hz  data  points  have  been  deleted  in  the  power-law  calculations 
in  Table  27.  In  general,  the  Biot  model  shows  slightly  better  agree¬ 
ment  with  the  data  over  the  frequency  range  than  does  Hamilton's  model. 
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TABLE  27 

POWER-LAW  FITS  TO  BOTTOM -ABSORPTION  COEFFICIENTS  (a.) 

A 

Panama  City 


Predicted 

aA 

Experimental  a. 

(dB/km) 

Leg 

(dB/km) 

Biot 

Hamilton 

C 

0.0177 

f1.81 

(Caswell,  1979) 

0.003  f1*86 

0.501  f1,0 

0.0028 

f1.75 

(Ingenito ,  1973) 

0.0290 

f1.59 

D 

0.0015 

f2.08 

(Caswell ,  1979) 

0.003  f1*86 

0.501  f1,0 

0.0090 

f  1. 71 

Note:  f  is  in  Hz 
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although  the  agreement  is  not  as  good  as  that  obtained  at  other  sites. 

As  may  be  seen  in  Table  27,  the  frequency  dependence  for  all  the  data 

(except  the  predictions  of  Hamilton)  is  closer  to  f  than  f  .  The 

difference  in  bottom  models  has  a  relatively  small  effect  on  the  value 

of  K,  and  n  as  may  be  seen  by  comparing  Caswell's  original  Leg  C 
A 

data  and  the  experimental  absorption  data  for  Leg  C  derived  here. 


CHAPTER  VII 


PREDICTION  OF  TRANSMISSION  LOSS 
SCOTIAN  SHELF  EXPERIMENTS 


7.1  Introduction 

The  ultimate  test  of  a  sediment  model  is  how  well  transmission 
loss  is  predicted  when  the  results  of  this  model  are  used  in  the 
propagation  model.  Thus,  in  this  chapter,  bottom  absorption-loss  (aA) 
is  predicted  based  on  the  best  available  sediment  data.  The  adiabatic 
normal-mode  propagation  model  is  then  run  using  the  bottom-absorption 
data  to  calculate  the  transmission  loss.  The  Scotian  Shelf  site  is 
interesting  for  this  purpose,  since  a  wide  variety  of  sediment  types 
and  bottom  structures  is  present.  In  addition,  excellent  sediment  data 
and  advice  as  to  its  interpretation  were  provided  by  the  Bedford 
Institute  of  Oceanography  of  Dartmouth,  Nova  Scotia. 

The  procedure  followed  is  similar  to  previous  chapters  in  that 
the  experiment  and  the  calculation  of  experimental  transmission  loss 
are  described  first.  The  prediction  of  transmission  loss  is  discussed 
next,  wherein  the  construction  of  a  geoacoustic  model  for  the  water- 
seabed  system  is  described.  Finally,  the  experimental  and  predicted 
loss  data  are  discussed  and  then  compared. 

7 . 2  Experiments 

During  September  of  1978,  experiments  were  performed  at  two 


sites  off  the  coast  of  Nova  Scotia.  Long-range  and  short-range 
propagation  and  seismic-refraction  experiments  were  performed  at 


134 


both  sites;  however,  only  the  long-range  experiment  at  one  site  will 
be  discussed  here.  Prior  to  the  September  experiment,  sub-bottom 
profiling  was  performed  along  the  propagation  track  using  the  Huntec 
Deep-Tow  System.  The  data  obtained  from  this  system  and  grab-sample 
data  obtained  from  the  Bedford  Institute  of  Oceanography  of  Dartmouth, 
Nova  Scotia,  provided  useful  information  describing  the  seabed  in  the 
test  area. 

At  the  test  site,  the  new  array  was  moored  in  water  57  m  deep. 

The  long-range  propagation  track  covered  a  sloping  bottom  with  the  water 
depth  reaching  approximately  265  m  near  the  maximum  range  of  125  km. 
Large  (456  g)  charges  were  dropped  at  nominal  depths  of  38  and  79  m  at 
range  intervals  of  1.6  km.  The  array  covered  approximately  the  bottom 
half  of  the  water  column  at  this  site.  The  bathymetry  and  sound- 
velocity  profiles  obtained  along  the  propagation  track  are  depicted  in 
Figure  30.  The  source  depths  are  indicated  by  the  clotted  lines  in  the 
figure. 


7 . 3  Calculation  of  Experimental  Transmission  Loss 

Signals  received  on  the  bottom  hydrophone  of  the  array  (1  m  above 
the  bottom)  were  analyzed  as  discussed  in  Chapter  III  to  determine  their 
energy  content  in  one-third-octave  bands  centered  at  frequencies  from 
25  to  800  Hz.  Source  spectrum  levels  for  the  456  g  shots  were  estab¬ 
lished  by  R.  C.  Hughes,  Defense  Research  Establishment  Atlantic, 
Dartmouth,  Nova  Scotia,  using  a  modification  of  the  Gaspin  and  Shuler 
(1971)  technique  [see  Hughes  (1976)]  to  obtain  improved  results  below 
the  bubble  frequency.  Hughes  verified  the  predicted  source  levels  for 
the  79  m  depth  by  comparing  them  with  actual  recorded  levels  for  3 
similar  charge  detonated  at  a  depth  of  77  m.  The  transmission  loss 
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Figure  30.  Depth  and  Sound  Velocity  Profiles — Scotian  Shelf 
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was  calculated  by  taking  the  difference  between  the  received  energy 
level  and  the  source  spectrum  level  in  the  desired  one-third-octave  band 
and  at  the  source-receiver  range  of  interest.  The  range  separation  was 
calculated  (see  Chapter  III)  using  the  shot  instant  time  (recorded  on 
the  shot  boat)  and  the  receive  time  to  obtain  the  travel  time;  the  range 
was  calculated  by  multiplying  the  travel  time  by  the  sound  speed  in  water . 

The  experimental  transmission  loss  measurements  for  the  test  site 
will  be  discussed  later  when  comparisons  with  the  predicted  losses  are 
made . 

7.4  Calculation  of  Predicted  Transmission  Loss 

7.4.1  Range-dependent  model.  Because  of  the  sloping  bottom 
found  at  the  site,  the  standard  normal-mode  propagation  model  does  not 
apply.  As  discussed  in  Chapter  II,  a  modification  to  this  model,  the 
adiabatic  normal-mode  model,  may  be  used  for  bottoms  having  a  moderate 
slope.  This  model  assumes  that  each  normal  mode  propagates  indepen¬ 
dently  with  no  exchange  of  energy  between  modes  (i.e.,  no  mode  coupling.) 
McDaniel  (1977)  has  shown  that  mode  coupling,  which  may  be  due  to 
surface  or  bottom  roughness,  inhomogeneities  within  the  water  column, 
or  severe  bottom  slopes,  increases  with  frequency.  Hence,  the  model 
will  have  definite  high-frequency  limitations. 

The  adiabatic  normal-mode  model  is  implemented  by  assuming  that 
the  propagation  track  may  be  divided  into  segments  for  which  the 
standard  normal-mode  theory  applies.  The  mode  excitation  may  be  found 
at  the  source  and  receive  positions  for  each  segment  and  then,  knowing 
the  mode  attenuation  along  the  segment,  the  transmission  loss  may  be 
calculated.  To  calculate  mode  attenuation,  estimates  of  the  attenuation 


due  to  scattering,  conversion  to  shear,  and  bottom  absorption  must  be 
obtained.  The  prediction  calculation  is  thus  the  reverse  of  that 
employed  in  Chapters  I  V  -  VI.  For  example,  bottom  absorption  will 
be  predicted  and  the  mode-attenuation  coefficient  calculated  from  the 
product  of  the  bottom-absorption  coefficient  and  the  energy  ratio. 

Thus,  both  the  water  column  and  ocean  bottom  must  be  well  characterized 
(i.e.,  a  geoacoustic  model  for  the  environment)  to  use  the  adiabatic 
normal-mode  model. 

7.4.2  Geoacoustic  model.  This  model  is  a  complete  environmental 
model  for  the  water-bottom  system  in  the  sense  that  all  the  parameters 
relating  to  acoustic  propagation  are  specified.  These  parameters 
include  the  sound  speed  in  the  water  and  its  variation  with  depth,  the 
compressional  and  shear  velocity,  the  density,  and  the  absorption 
coefficients  for  compressional  and  shear  waves  as  a  function  of  depth 
in  the  bottom.  As  noted  previously,  the  bottom  is  assumed  to  consist 
of  discrete  layers  so  that  the  layer  thickness  and  the  other  parameters 
must  be  specified  for  each  bottom  layer. 

Because  of  the  Huntec  profiling  and  the  grab-sample  data 
available  for  the  Scotian  Shelf,  the  formulation  of  a  well-documented 
geoacoustic  model  was  possible.  In  addition,  seismic-refraction  and 
dispersion  analysis  provided  confirmation  (and,  in  some  cases,  the  only 
data)  for  the  chosen  model. 

The  Huntec  system  consists  of  an  electrodynamic  source  towed  at 
mid-depth  and  the  associated  shipboard  signal  processing  equipment. 

It  has  been  shown  by  Parrott  et  al.  (1979)  that  the  system  provides  a 
reliable  indication  of  surficial  sediment  type  due  to  a  statistical 
analysis  of  the  energy  returned  from  the  bottom.  The  presence  and 
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identity  of  deeper  (up  to  50  m  or  so)  sub-bottom  reflectors  is  readily 
apparent  on  the  system  output.  In  many  cases,  surficial,  intermediate, 
and  bedrock,  layers  were  identified  from  the  records.  Where  grab-samples 
were  available  near  the  propagation  track,  the  mean  grain-size  obtained 
from  this  data  was  used  in  the  sediment  models  to  predict  bottom 
absorption.  In  other  cases,  the  Huntec  records  were  used  as  an 
indication  of  sediment  mean  grain-size.  The  findings  of  King  (1970) 
provided  supporting  data  on  the  sediment  types  at  the  site. 

Sound  speeds  were  established  using  a  regression  equation 
derived  by  Anderson  (1974)  which  used  mean  grain-size  as  the  independent 
variable.  The  velocity  c (m/sec)  at  a  temperature  of  20°C  and  a 
pressure  of  one  atmosphere  was  calculated  from  the  expression: 

c  =  -68.156  M  +  3.053  M  2  +  1874.4 

z  z 

These  velocities  were  then  corrected  to  the  bottom  temperature  and 
depth. 

The  intermediate  layers  were  identified  by  the  appearance  of  the 
fine  structure  on  the  Huntec  records.  A  mean  grain-size  was  assigned 
based  on  this  identification  and  the  absorption  was  calculated  using 
the  sediment  models.  The  work  of  Hamilton  (1975)  was  followed  in 
determining  a  depth  correction  for  the  calculated  absorption  values. 

In  those  cases  where  the  bedrock  layer  could  not  be  identified 
from  the  Huntec  record,  estimates  of  the  depth  to  this  layer  and  the 
probable  type  were  obtained  from  King  and  MacLean  (1976).  Tertiary 
and  cretaceous  formations  were  identified  as  the  bedrock  structure 
along  the  track.  King  estimated  that  compressional  sound  speeds  of 


2000  to  2300  m/sec  would  exist  for  these  structures;  hence,  the 
bedrock  is  classified  as  more  of  a  semi-consolidated  sediment. 

Based  on  the  preceding  considerations,  the  bottom  structure  at 
the  site  was  modeled  as  shown  in  Figure  31.  For  use  in  the  adiabatic 
normal-mode  model,  the  bottom  in  the  figure  was  broken  down  into  the 
segments  identified,  by  range  from  the  array,  in  Table  28.  At  the  site 
and  along  the  track,  sediment  identification  was  based  on  grab-sample 
data  as  noted  in  the  table.  Beneath  the  array,  the  sediments  consisted 
of  gravel  and  sand.  Because  the  regression  equation  for  calculating 
the  velocities  does  not  apply  to  such  coarse  sediments,  seismic- 
refraction  and  dispersion  analyses  were  used  to  establish  the  first 
layer  velocity  for  the  zero  and  11  km  segments.  The  19  km  segment 
consisted  of  slightly  finer  sediments;  hence,  the  velocity  was  reduced 
slightly  from  that  of  the  previous  segments.  The  portion  of  the  track 
from  27  through  95  km  was  modeled  using  velocities  and  absorption 
coefficients  predicted  by  the  regression  equation  and  the  Biot  sediment 
model  using  mean  grain-size  data.  At  the  opposite  end  of  the  track 
(117  km) ,  the  Huntec  records  showed  that  the  bottom  consisted  of  a  thin 
clay  layer  over  a  pocket  of  silt  and  diffused  gas.  A  mean  grain-size 
was  assumed  from  the  Huntec  data  for  the  clay.  The  absorption  and 
velocity  for  the  gassy  silt  were  obtained  by  modifying  the  fluid 
parameters  in  the  Biot  model  according  to  an  assumed  concentration  of 
methane  gas  in  seawater.  With  this  approach,  the  pore  fluid  was 
assumed  to  be  a  suspension  of  small  methane  bubbles  in  seawater  so 
that  the  equation  of  Wood  (1944)  applied. 

That  is,  the  density  and  compressibility  of  the  mixture  were 


obtained  from  the  volume  fraction  of  methane  and  seawater  present. 


SEABED  PARAMETERS  -  SCOTIAN  SHELF 
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TABLE  28  (continued) 
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Values  were  obtained  from  Table  29  for  0.1  and  0.01%  concentrations  of  methane. 
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Thus,  if  V  is  the  fraction  of  methane  present,  the  density  (P  ), 
m  mv 

bulk  modulus  (k  ),  and  sound  velocity  are  obtained  as  follows: 
mw 


P 

mw 


VP  +  (1  -  V  )P 

mm  m  w 


k 

mw 


k  k 
m  w 


V  k 
m  w 


+  (1  - 


V  )k 


m 


indicate  quantities  of  methane  and 

water,  respectively.  The  effect  of  changing  the  fraction  of  methane 

in  the  mixture  may  be  seen  in  Table  29.  The  values  for  K  ,  n  and 

A 

C  were  obtained  by  using  p  and  k  as  the  fluid  density  and 

mw  mw 

bulk  modulus  in  the  Biot  sediment  model.  As  will  be  discussed  later, 
adiabatic  normal-mode  runs  were  made  using  the  Biot  model  parameters 
from  Table  29  for  concentrations  of  0.1  and  0.01  percent  methane  within 
the  segment  at  117  km. 


and 


A- 

c  «  /-SK 
mw  P 

mw 


where  the  subscripts  m  and  w 


7.5  Comparison  of  Measured  and  Predicted  Transmission  Loss 

The  experimental  and  predicted  transmission  loss  data  are  plotted 
in  Figures  32  (38  ra  source  depth)  and  33  (79  m  source  depth)  for 
frequencies  of  25,  80,  250,  and  800  Hz.  The  discussion  to  follow  begins 
with  a  paragraph  relating  the  experimental  data  to  environmental 
conditions;  the  predicted  loss  is  discussed  next  in  terms  of  the 
geoacoustic  model,  and  finally,  possible  reasons  for  agreement  and 
disagreement  between  the  measured  and  predicted  loss  are  considered. 


EFFECT  OF  METHANE  CONCENTRATION 


TRANSMISSION  LOSS  (dB) 


TRANSMISSION  LOSS  (dB) 


Figure  33.  Comparison  of  Experimental  and  Predicted 

Transmission  Loss-Scotian  Slielf,  79  m  Depth 
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7.5.1  Experimental  transmission  loss.  Propagation  at  the 
Scotian  Shelf  site  was  largely  controlled  by  the  sound-velocity  profiles 
which  existed  there.  The  loss  measured  for  the  shallow  source  (Figure 
32)  reflects  this  best,  as  charges  were  dropped  over  all  but  the  first 

2  km  of  the  track.  In  Figure  30,  it  can  be  seen  that  out  to  a  range  of 
20  km,  the  shot  is  detonated  in  a  region  where,  due  to  the  decrease  of 
sound  velocity  with  depth  (a  negative  profile) ,  the  sound  rays  will  be 
refracted  into  the  bottom.  This  resulted  in  bottom-dominated  propaga¬ 
tion  in  this  region  and  accounted  for  the  large  increase  in  loss  as  the 
range  increased. 

Beyond  20  km,  the  sound-velocity  profiles  in  Figure  30  indicate 
the  presence  of  a  sound-velocity  minimum.  This  type  of  profile  will 
cause  the  energy,  for  a  source  detonated  near  the  axis  of  the  minimum, 
to  be  channeled  with  little  interaction  with  the  surface  or  bottom. 
Losses  in  the  channel  are  expected  to  be  near  cylindrical  spreading  in 
magnitude.  The  losses  were,  in  fact,  near  cylindrical  spreading  for 
25  and  80  Hz  beyond  20  km,  as  both  source  depths  were  near  the  channel 
axis  (note  the  dotted  lines  in  Figure  30) .  A  last  point  of  interest  is 
near  117  km  where  the  pocket  of  silt  with  diffuse  gas  is  indicated  in 
Figure  31.  An  increase  in  loss  is  shown  in  Figure  32  in  this  same 
region,  particularly  at  25  Hz,  and  In  Figure  33  for  all  frequencies  to 
some  extent.  In  general,  the  loss  appeared  to  increase,  for  both  source 
depths,  as  the  frequency  increased  from  25  to  800  Hz. 

7.5.2  Predicted  transmission  loss.  Using  the  sound-velocity 
profiles  from  Figure  30  along  with  the  seabed  parameters  of  Table  28, 
the  adiabatic  normal-mode  computations  of  transmission  loss  were 
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performed.  As  noted  in  the  table,  the  bottom  is  specified  by  12 
segments.  Over  the  first  20  km  of  the  propagation  track  in  Figure  32, 
the  bottom  is  well  characterized  by  the  first  three  segments,  and  high 
losses  are  predicted  for  the  coarse  sediments  in  this  region.  In  the 
region  between  20  and  95  km,  the  finer  silt  and  sand  sediments  modeled 
in  Table  28  cause  the  rate  of  loss  to  be  somewhat  reduced  from  that 
encountered  over  the  early  portion  of  the  track.  The  erratic  behavior 
shown  at  800  Hz  reflects  the  extreme  sensitivity  to  changing  environ¬ 
mental  parameters  encountered  as  the  frequency  increases.  The  segment 
at  117  km  reflects  profoundly  the  effect  of  the  silt-gas  sediments. 
Methane  concentrations  of  0.1  and  0.01  percent  are  used  since  the  true 
concentration  is  unknown.  The  effect  of  this  gas  is  to  reduce  the 
sound  velocity  and  increase  the  attenuation  within  the  sediment. 

Because  the  adiabatic  model  uses  interpolation  between  segments,  the 
effects  of  the  gassy  sediment  are  indicated  by  increased  transmission 
loss  from  95  km  to  the  end  of  the  track.  The  main  cause  of  the 
increased  loss  is  the  low  sound  velocity.  A  low-velocity  layer  is 
formed, as  in  the  Corpus  Christi  discussion,  which  traps  energy  in  the 
layer  and  reduces  the  pressure  field  in  the  water  column. 

The  transmission  loss  predicted  using  the  Hamilton  method  for 
sediment  absorption  is  also  shown  in  Figure  32. 

7.5.3  Comparison.  In  both  figures,  the  agreement  is  excellent 
for  frequencies  of  25,  80,  and  250  Hz.  At  25  Hz,  the  propagation  model 
(Biot  sediment  model)  accurately  predicts  the  rapid  increase  in  loss 
occurring  at  ranges  less  than  20  km.  As  discussed  earlier,  the  bottom 
absorption  dominates  this  loss  since  most  of  the  energy  is  refracted 


149 


into  Che  bottom.  The  reduced  loss,  occurring  at  ranges  greater  than 
20  km  due  to  the  sound  channel,  is  also  predicted  with  reasonable 
accuracy  out  to  a  range  of  90  km  or  so.  At  this  range,  the  experimental 
transmission  loss  tends  to  increase  somewhat  due  to  the  basin  where  the 
gassy  sediments  are  found.  The  loss  predicted  for  the  two  assumed 
methane  concentrations  appears  to  bracket  the  observed  variation  in 
transmission  loss.  That  is,  a  reduction  of  the  second  bottom  layer 
sound  velocity  to  942  m/sec  (0.1  percent  concentration)  in  the  model 
causes  an  excess  of  energy  to  be  trapped  in  the  layer.  The  higher 
velocity  (1411  m/sec)  for  the  0.01  percent  concentration  has  very  little 
effect.  Possibly,  an  intermediate  concentration  or  a  more  detailed 
model  (in  range  segments)  would  give  better  results. 

The  loss  predicted  at  80  Hz  gives  excellent  agreement  over  all 
of  the  track  if  the  methane  concentration  is  correctly  assumed  to  be 
0.01  percent.  The  experimental  data  shows  little  of  the  effects  (for 
either  source  depth)  at  117  km  that  were  noted  at  25  Hz,  so  the  methane 
could  be  ignored  with  probably  the  same  result. 

As  the  frequency  increases  to  250  Hz,  the  rapi-'  increase  in  loss 
over  the  first  20  km  of  the  track  is  still  predicted.  At  ranges 
greater  than  20  km,  the  agreement  for  the  shallov  source  is  not  quite 
as  good  as  for  the  deep  source.  In  the  region  around  117  km,  agreement 
is  still  good  for  both  source  depths,  assuming  a  0.01  percent  concen¬ 
tration.  The  shallow  source  and  the  reduced  bottom  penetration  (due  to 
the  higher  frequency)  appear  to  result  In  less  trapping  of  energy  in 
the  low-velocity  layer  for  the  0.1  percent  concentration.  The 
deep  source  still  predicts  losses  in  excess  of  100  dB  in  the  gassy 
region  for  the  same  concentration. 


The  predicted  losses  at  800  Hz  are  in  poor  agreement  with  the 
measured  data.  While  the  frequency  dependence  of  the  bottom  absorption 
appears  to  be  satisfactory  over  the  lower  frequencies,  an  excess  of 
attenuation  is  predicted  at  800  Hz.  The  only  region  with  good  agreement 
is  the  section  between  40  and  90  km  for  the  shallow  source.  Extreme 
effects  are  noted  beyond  90  km  for  both  concentrations.  Since  this 
behavior  is  not  consistent  with  the  low  frequency  performance,  something 
other  than  the  bottom  model  is  assumed  responsible.  It  appears  that 
these  predicted  variations  are  due  to  the  changing  temperature  structure 
shown  in  Figure  30,  and  particularly,  the  depth  of  the  isovelocity 
layer.  In  addition,  at  this  frequency,  mode  coupling  can  no  longer 
be  ignored  and  its  inclusion  in  the  model  will  have  some  effect 
(although  unknown  at  this  time).  Thus,  the  data  at  800  Hz  does  not 
represent  a  valid  test  of  our  ability  to  predict  bottom  absorption  as 
the  propagation  model  itself  is  not  adequate  at  these  frequencies  and 
water  depths. 

In  Figure  32,  transmission  loss  is  also  predicted  using 

Hamilton's  method  for  calculating  sediment  absorption.  The  mean 

grain-sizes  in  Table  28  were  used  in  Hamilton's  equations  to  obtain 

a,  ,  and  these  values  were  then  used  in  the  adiabatic  normal-mode 
A 

propagation  model.  At  117  km,  the  gassy  silt  was  modeled  with  the 
Biot  values  (0.01  percent  concentration)  for  absorption  and  velocity, 
although  the  upper  clay  layer  was  modeled  using  Hamilton's  equation. 

As  can  be  seen  in  Figure  32,  the  transmission  loss  predicted  at  23  Hz 
is  much  greater  than  that  measured  for  ranges  beyond  15  km  or  so.  At 
80  Hz,  the  predictions  are  only  4  to  5  dB  low,  and  at  250  Hz,  the 
situation  reverses  with  predictions  somewhat  less  than  that  measured. 
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The  loss  predicted  at  800  Hz  parallels  the  Biot  results  but  is  4  to 
5  dB  greater.  Overall,  the  results  using  Hamilton's  method  give 
poorer  agreement  than  those  using  Biot's  theory. 


CHAPTER  VIII 


SUMMARY  AND  CONCLUSIONS 

The  degree  to  which  transmission  loss  may  be  accurately  predicted 
is  dependent  upon  both  a  detailed  knowledge  of  the  ocean  environment  and 
the  use  of  propagation  models  which  faithfully  represent  the  physical 
processes  at  work.  The  purpose  of  this  study  has  been  to  show  that, 
given  a  comprehensive  propagation  model  (the  normal-mode  model) , 
transmission  loss  may  be  predicted  using  an  environmental  model  whose 
bottom  effects  are  predicated  on  two  sediment  models.  As  a  test  of  the 
accuracy  of  the  theories,  experimental  bottom  absorption,  calculated 
from  mode-attenuation  coefficients  using  a  normal-mode  model,  has  been 
compared  with  the  absorption  predicted  by  the  Biot  sediment  model. 

This  comparison  was  also  performed  for  sediment  absorption  calculated 
using  the  procedures  developed  by  Hamilton. 

Values  of  bottom  absorption  obtained  from  experimental  data 
acquired  at  several  test  sites  were  compared  with  the  bottom  absorption 
predicted  by  the  Biot  and  Hamilton  models.  The  experimental  bottom 
absorption  was  calculated  through  the  relation  aA  =  6/y  ,  where 

the  mode-attenuation  coefficient  (5)  was  obtained  experimentally.  The 
energy  distribution  (y)  was  obtained  by  calculations  performed  using 
the  normal-mode  program.  The  accuracy  of  this  calculation  for  a. 


rested  upon  the  validity  of  the  bottom  model  developed  to  represent 
the  test  site.  A  key  component  of  the  discussion  within  these  chapters 
was  the  variation  of  the  experimental  and  predicted  bottom  absorption 
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with  frequency.  This  frequency  dependence  constitutes  one  of  the 
basic  differences  between  the  Biot  and  Hamilton  predictions  of  sediment 
absorption. 

The  Scotian  Shelf  data  were  used  for  testing  the  transmission 
loss  predictions  using  the  range-dependent  propagation  model  and  the 
two  sediment  models.  A  wide  variation  in  sediment  type  from  gravel  to 
clay  was  present  along  the  propagation  track  and  excellent  data 
characterizing  the  seabed  was  available. 

8.1  Discussion  of  Results 

The  calculation  of  experimental  bottom  absorption  was  based  on 
mode-attenuation  data,  where  sediment  absorption  was  the  dominant  loss 
mechanism,  and  an  environmental  model  which  accurately  represented  the 
water-bottom  system.  The  mode-attenuation  data  have  been  discussed  in 
some  detail  in  the  relevant  chapters,  as  have  the  environmental  data. 
However,  some  additional  discussion  of  the  environment  is  appropriate. 
In  addition,  the  sediment  models  and  the  parameters  required  to  use 
them  are  discussed  and  compared. 

8.1.1  Environmental  effects.  The  calculated  bottom  absorption 
is  dependent  upon  the  environmental  model  through  the  energy  distribu¬ 
tion  parameter  7  .  That  is,  any  change  in  the  sound-velocity  profile 
within  the  water  column  or  within  the  bottom  will  affect  the  mode 
amplitude  distribution  in  the  water  and  bottom  and,  hence,  the  energy 
distribution.  These  effects  are  frequency  dependent  in  that  the  high 
frequencies  are  influenced  much  more  than  the  low  frequencies.  For 
any  sound  velocity  profile,  most  of  the  energy  will  be  concentrated  in 
the  water  column  at  high  frequencies  (750  Hz  or  greater).  Of  the 
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remaining  energy  (that  in  the  bottom) ,  99  percent  is  found  within  the 
uppermost  bottom  layer.  This  is  in  contrast  to  the  situation  at  low 
frequencies  (50  Hz  and  less)  where  the  energy  in  the  uppermost  layer 
may  be  only  40  percent  or  so  of  the  total  bottom  energy.  Because  of 
this,  a  small  change  in  the  sound-velocity  gradient  (e.g.,  from  a  weak 
sound  channel  to  a  positive  gradient)  can  cause  a  decrease  in  y^^ 
and,  consequently,  a  significant  increase  (50  percent  or  so)  in  the 
absorption  at  high  frequency;  essentially  none  occurs  at  low 
frequency.  As  the  sound  waves  propagate  along  a  track  where  such 
changes  occur,  only  net  effects  are  observed.  Hence,  by  choosing  some 
mean  profile,  these  effects  can  be  modeled.  The  situation  at  Daytona 
Beach — Site  1  represents  this  type  of  problem,  and  by  choosing  a 
representative  sound-velocity  profile,  the  changing  sound  speed  profile 
effects  are  minimized.  In  all  the  other  test  data,  except  the 
range-dependent  case  on  the  Scotian  Shelf,  a  near-isovelocity  sound 
speed  profile  was  present. 

The  bottom  parameters  affecting  the  experimental  bottom 
absorption  are,  for  the  three-layer  model,  the  velocity  (c^)  and 
thickness  (H^)  of  the  first  bottom  layer,  and  the  velocity  (c^)  of 
the  second  bottom  layer.  The  effect  of  certain  of  these  parameters 
has  been  noted  with  regard  to  the  Panama  City  data.  The  general 
behavior  of  these  parameters  may  be  inferred  through  a  more  detailed 
analysis  of  the  Panama  City — Leg  C  bottom  absorption  data.  The 
bottom  model  used  for  this  site  is  listed  (model  number  4)  in 
Table  30  (a)  along  with  four  other  possible  bottom  models  and  the 
consequent  bottom  absorption  (a.)  for  frequencies  of  50,  400,  and 
750  Hz.  Sensitivity  parameters  SC^  ,  SH^  ,  and  SC2  are  calculated 
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in  Table  30  (b)  by  dividing  the  change  in  by  the  percentage 

change  in  the  parameter  of  interest,  the  other  parameters  remaining 
constant.  The  two  models  involved  for  each  comparison  (e.g.,  two  and 
three  for  SC^  )  and  the  resulting  sensitivity  to  the  three  parameters 
at  each  frequency  are  noted.  As  can  be  seen,  at  50  Hz,  the  absorption 
is  somewhat  more  sensitive  to  than  the  other  parameters.  At  the 

mid-frequency  (400  Hz)  and  the  highest  frequency  (750  Hz),  the 
absorption  is  extremely  sensitive  to  .  An  analysis  of  the 
sensitivity  of  and  n  would  show  a  similar  relation  (i.e.,  most 

sensitive  to  c^  )  as  has  a  large  effect  over  a  great  proportion 

of  the  frequency  range.  The  negative  values  for  SH^  indicate  that, 
as  the  layer  thickness  decreases,  the  absorption  increases.  The  value 
selected  for  c^  in  the  bottom  model  for  each  site  is  based  on  the 
best  available  information;  for  the  Daytona  Beach  and  Corpus  Christi 
sites,  c^  was  determined  by  dispersion  analysis,  and  for  the  Panama 
City  site,  c^  was  based  on  an  in-situ  measurement  with  the  coring 
apparatus.  This  value  will  be  most  reliable  for  short-range 
measurements. 

8.1.2  Predicted  bottom  absorption.  The  ability  to  predict 
bottom  absorption  is  dependent  upon  the  adequacy  of  the  sediment  data 
and  the  validity  of  the  sediment  model.  However,  even  the  best  of 
sediment  information  will  describe  only  the  upper  few  meters  of  the 
ocean  bottom  and  waves  propagating  at  the  lower  frequencies  will 
frequently  penetrate  much  deeper  than  this. 

The  ideal  type  of  sediment  data  consists  of  long  cores,  taken 
every  kilometer  or  so  along  the  shot  tracks.  This  enables  a  complete 
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description  of  the  sediments  and  their  variation  with  range  and  depth 
to  be  made.  As  noted  in  the  respective  analysis  chapters,  this  type 
of  data  was  generally  not  available.  Of  the  sites  discussed,  the 
short-range  Corpus  Christi  site  has  the  best  sampling  of  cores;  two 
cores  were  obtained  near  the  site,  and  two  more  were  spaced  along  the 
short-range  track.  In  addition  to  the  standard  analysis  of  these  cores, 
the  sound  velocity  through  the  core  was  also  measured.  The  combination 
of  Huntec  profiling  and  ample  grab-samples  resulted  in  good  sediment 
data  for  the  Scotian  Shelf  site.  The  Panama  City  and  Daytona  Beach 
sediments  were  not  characterized  as  well.  Adequate  coring  data  were 
available  at  the  deeper  site  (31  m)  off  Panama  City,  but  none  was 
available  along  the  track  at  the  deep  site  nor  for  the  shallow  site 
and  shot  track.  The  shallow  site  (Site  1)  off  Daytona  Beach  was 
sampled  with  cores  and  grab-samples  around  the  periphery  of  the 
short-range  legs,  while  the  deeper  site  was  characterized  by  grab- 
samples  at  either  end  of  the  long-range  shot  track.  In  most  cases, 
except  for  the  Scotian  Shelf,  the  sediment  parameters  assumed  are 
most  representative  of  the  short-range  propagation  tracks.  These 
parameters  (  Mz  and  0^  )  are  listed  in  the  table  of  summary  parameters. 
Table  31,  along  with  other  parameters  to  be  discussed  later. 

While  the  best  available  data  were  used  as  inputs  for  the 
sediment  models,  all  of  the  information  desired  was  not  available  at 
every  site.  Thus,  the  validity  of  the  sediment  model  is  not  so  much 
a  question  of  its  credibility  as  a  question  of  the  validity  of  its 
input  and  peripheral  calculations.  Both  sediment  models  have  achieved 
credibility  as  discussed  in  Chapter  II.  Thus,  the  critical  input 
parameters  will  be  discussed  and  then  certain  effects  embodied  in  the 


158 


159 


Biot  theory  will  be  considered.  These  effects  are  the  frequency 
dependent  viscous  losses  and  the  frame  losses. 

The  most  sensitive  of  the  input  parameters  is  the  permeability 
which  is  calculated,  Equation  (46),  from  the  mean  grain  diameter  (in 
mm)  and  its  standard  deviation.  The  permeability  is  smallest  for  the 
fine  grain  sizes,  and  poorer  sorting  (large  a )  reduces  the  permea¬ 
bility.  Consequently,  a  poorly  sorted  sediment  has  the  absorption 
characteristics  of  a  finer  sediment.  Hovem  (1980)  has  derived  an 
absorption  relation  using  a  modified  suspension  theory  in  order  to 
include  the  effects  of  sorting.  As  noted  in  Chapter  II,  sorting  is 
handled  within  the  permeability  calculation  in  an  empirical  manner 
and  has  the  greatest  effect  on  the  coarser  grain-sizes. 

Within  the  Biot  theory,  losses  due  to  the  inelastic  nature  of 

the  granular  frame  are  included  by  making  the  shear  and  bulk  modulus 

(for  the  frame)  complex  functions.  This  is  expressed  by  the  pair  of 

equations  for  kfa  and  y  [Equation  (39)].  These  are  calculated,  as 

discussed  by  Stoll,  from  experimental  values  for  the  log  decrement 

(Ar  and  A  )  and  the  longitudinal  and  shear  velocity  for  the  dry, 

granular  frame  material  (v^,  and  vg).  These  steps  involve  Equations 

(40)  through  (44)  and,  based  on  Stoll's  work,  require  different  values 

for  A„  ,  A  ,  and  a  ,  depending  on  the  grain-size.  Thus,  to  provide 
c.  s 

consistency,  values  of  0.15,  0.20,  and  1.2  were  assumed  for  A^  ,  Ag  , 
and  ol  ,  respectively,  for  mean  grain-sizes  less  than  2.5  phi,  and 
0.45,  0.60,  and  3.0  for  greater  than  2.5  phi.  Above  4  phi,  the 

main  effect  of  selecting  the  larger  values  is  to  increase  the  magnitude 
of  the  attenuation  and  not  the  frequency  dependence.  Stoll  reports 


good  agreement  with  other  fine  sediment  data,  but  the  data  reported 


here,  although  in  good  agreement,  is  not  considered  extensive  enough 
to  provide  unquestioned  verification  of  Stoll's  approach  to  this 
aspect  of  the  theory. 


The  effect  of  losses  due  to  frame  inelasticity  (e.g.,  inter¬ 
granular  friction)  may  be  evaluated  by  removing  this  mechanism, 
setting  A^  =  Ag  =  0  ,  and  comparing  the  absorption  with  that  predicted 
using  the  values  noted  for  the  log  decrement.  This  was  done  for  mean 
grain-sizes  ranging  from  0  to  7  phi  in  1  phi  increments  (o  =  0.5)  over 
a  frequency  range  of  50  to  750  Hz.  The  results  indicate  that  frame 
losses  play  a  dominant  role  for  ^  2  phi  ,  based  on  an  arbitrary 
threshold  defined  by  absorption  values  within  10  percent  of  the  normal 
case.  This  result  is  somewhat  frequency  dependent  as  the  equality  in 
the  above  relation  holds  for  frequencies  of  100  Hz  or  less.  As 
expected,  this  is  in  agreement  with  Stoll. 

The  contribution  from  viscous,  flow-related  losses  was  evaluated 
over  the  same  range  of  grain-size  and  frequency  by  forcing  the  function 
F  ,  in  Equations  (30)  and  (31),  to  be  unity.  The  results  show  that 
this  type  of  loss  dominates  for  <  2  phi  ,  based  on  the  same  criteria 
as  above.  Some  frequency  dependence  was  noted  for  =  1  ,  as  the  flow 
losses  dominated  at  400  Hz  and  above. 


8.1.3  Accuracy  of  predicted  absorption.  In  Chapters  IV  through 
VI,  absorption  was  predicted  using  the  methods  of  Biot  and  Hamilton. 

The  results  were  compared  with  the  experimental  absorption  calculated 
from  mode-attenuation  coefficients.  These  comparisons  were,  in  general, 
good,  but  in  recalling  this  data  (Figures  16,  21,  27,  and  29),  it  is 
obvious  that  the  predicted  absorption  agrees  better  with  the  experi¬ 
mental  results  for  some  sites  than  for  others.  It  is  also  apparent 


that  the  Biot  model  appears  to  give  better  results  than  the  method  of 
Hamilton.  However,  if  some  numerical  measure  of  closeness  could  be 
calculated,  the  comparisons  would  have  more  meaning.  This  was  done  by 
calculating  the  root-mean-square  (RMS)  difference  over  the  applicable 
frequency  range  between  the  experimental  absorption  values  and  the 
values  predicted  by  the  models.  Thus,  from  Tables  4,  11,  19,  and  25, 
experimental  values  of  a were  obtained,  and  using  the  power  law  form 
of  the  predicted  absorption  from  Table  31,  differences  were  calculated 
for  each  frequency.  The  RMS  values  of  these  differences  were  then 
calculated  for  each  set  of  data  and  are  listed  in  Table  31. 

The  magnitude  of  the  RMS  difference  is  a  measure  of  how  well  the 
predicted  absorption  fits  the  experimental  data.  Thus,  in  Table  31, 
the  site  with  the  best  agreement  is  seen  to  be  the  Corpus  Christi  site, 
and  the  site  with  the  worse  agreement  is  Leg  C  at  Panama  City.  The 
data  has  been  edited  by  deleting  some  high-frequency  data  as  discussed 
previously  so  that  for  the  data  considered,  the  dominant  loss  mechanism 
is  volume  absorption  in  the  sediments.  However,  the  poor  agreement  at 
the  higher  frequencies  is  still  the  greatest  contributor  to  the  large 
RMS  difference  for  Leg  C.  In  comparing  the  RMS  values  obtained  using 
the  Biot  and  Hamilton  approaches,  it  is  obvious  that,  in  all  cases,  the 
Biot  model  predictions  give  the  best  fit.  It  is  of  interest  to  note 
that  the  Corpus  Christi  site,  at  which  the  best  agreement  between 
measured  and  predicted  absorption  was  obtained,  is  also  the  site  at 
which  the  best  bottom  data,  in  the  form  of  cores  along  the  shot  track, 
was  obtained. 

The  frequency  dependence  of  bottom  absorption  is  another  point 
of  interest  illustrated  in  Table  31.  The  experimental  data  collected 
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by  Hamilton  over  a  wide  range  of  frequencies  show  a  first  power 
dependence  on  frequency;  however,  as  has  been  discussed,  this  point 
is  disputed  by  many  researchers.  As  can  be  seen  in  Table  31,  power 
law  fits  to  the  experimental  data  show  a  first  power  dependence  only 
for  the  long-range  shot  track  at  Corpus  Christi.  Some  variation  does 
exist,  but  in  all  other  cases,  the  frequency  dependence  of  the 
experimental  data  is  much  greater  than  the  first  power.  The  predictions 
of  the  Biot  model  tend  to  show  either  unity  or  something  approaching 
a  second  power  frequency  dependence. 

8 . 2  Conclusions 

As  noted  previously,  the  ultimate  goal  of  developing  a  valid 
sediment  model  is  to  improve  transmission  loss  predictions.  Sediment 
absorption  is  one  loss  mechanism  among  several  which  must  be  modeled 
effectively  if  transmission  loss  predictions  are  to  be  improved.  A 
site  off  the  coast  of  Nova  Scotia  was  selected  to  evaluate  the  Biot 
and  Hamilton  models  for  computing  sediment  absorption  within  the 
adiabatic  normal-mode  propagation  model.  The  Nova  Scotia  site  presented 
difficulties  due  to  its  range  dependence,  but  it  was  also  very  well 
characterized  by  a  number  of  data  sources.  This  data  implies  that 
sediment  absorption  was  the  dominant  loss  mechanism  along  the 
propagation  track.  The  results  of  the  transmission  loss  predictions 
for  this  site  in  Chapter  VII  appear  to  confirm  this,  as  good  agreement 
was  obtained  for  frequencies  from  25  to  250  Hz.  The  results  at  800  Hz 
were  much  poorer  for  reasons  discussed  in  Section  7.5.3. 
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The  processes  involved  in  constructing  the  seabed  model  for 
Che  propagation  track  and  the  subsequent  results  lead  to  several 
conclusions,  some  of  which  relate  to  the  sediment  model,  while  others 
reflect  upon  the  propagation  model  itself. 

(1)  The  Biot  model,  with  appropriate  inputs,  provides  good 
estimates  of  the  bottom  absorption  needed  to  model 
propagation  loss.  The  Hamilton  model  performs  well 
for  fine-grained  sediments,  but  does  not  perform 
adequately  for  coarse-grained  sediments  due  to  its 
linear  dependence  upon  frequency. 

(2)  It  is  extremely  important  to  know  as  much  about  the 
deeper  sediment  layers  as  about  the  surficial  layer. 

The  depth  of  interest  is  approximately  one  wavelength 
so  that,  if  the  lowest  frequency  considered  is  25  Hz 

as  for  the  Nova  Scotia  site,  then  the  bottom  parameters 
must  be  specified  to  a  depth  of  60  m  or  so. 

(3)  The  combination  of  high  resolution  sub-bottom  profiling 
and  deep  cores  will  provide  this  type  of  data.  For 
the  Nova  Scotia  site,  grab-samples  rather  than  cores 
were  available;  however,  the  site  was  one  which  had 
been  extensively  studied  by  the  Bedford  Institute  of 
Oceanography.  The  identification  of  a  sediment  type 
from  Huntec  data,  or  actual  grab-sample  data,  gave  a 
mean  grain-size  which  was  then  used  in  the  Biot  sediment 
model  to  obtain  the  absorption  values  for  propagation 
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(4)  A  valuable  attribute  of  the  Biot  model  is  the  ability 
to  modify  both  fluid  and  frame  parameters.  This  was 
of  particular  value  when  modeling  the  sediments  con¬ 
sisting  of  silt  containing  diffuse  gas.  At  other 
sites,  the  bulk  modulus  for  the  frame  material  was 
modified  according  to  its  mineralogy  (i.e.  ,  the 
calcium  carbonate  content) . 

(5)  Frame  losses  require  further  studies  to  verify  Stoll's 
use  of  the  log  decrements  in  establishing  complex  bulk 
and  shear  moduli  for  fine-  and  coarse-grained  sediments. 

(6)  Modeling  of  sound  propagation  through  gas-filled 
sediments  must  be  improved.  These  sediments  are  found 
in  both  deep-sea  and  shallow-water  areas.  Mclver  (1974) 
has  summarized  tests  for  gas  content  performed  on  35 
cores  obtained  through  the  Deep  Sea  Drilling  Project 
(DSDP)  and  found  that  the  gas  is  98  to  100  percent 
methane  in  concentrations  (for  individual  cores) 
ranging  from  0.11  to  21.7  percent  by  volume.  The 
cores  were  obtained  on  the  Oregon  Shelf,  in  the 
Aleutian  Trench,  the  Bering  Sea,  and  the  Red  Sea. 

By  averaging  the  concentrations  obtained  on  cores  from 
the  same  region,  the  concentrations  ranged  from  1.0  to 
15.7  percent.  Thus,  the  concentrations  assumed  for  the 
Scotian  Shelf  site  are  possibly  too  low. 

(7)  Propagation  modeling  at  800  Hz  for  the  Scotian  Shelf 
and  other  areas  with  range-dependent  environments 
necessitates  the  use  of  a  model  which  includes  mode 
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coupling.  To  take  advantage  of  this  added  sophistication, 
greater  detail  in  the  environmental  model  may  be 
required. 

(8)  The  more  general  case  of  propagation  modeling  requires 
the  inclusion  of  other  loss  mechanisms.  For  example, 
sites  having  near-surface,  high-velocity  bedrock  will 
require  a  solid-bottom  model  and  absorption  coefficients 
for  both  compressional  and  shear  waves.  Very  rough 
bottoms  or  high  seas  will  require  that  scattering  losses 
and/or  mode  coupling  be  included  at  even  lower  frequencies. 
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APPENDIX  A 


NORMAL-MODE  INTEGRAL  EVALUATION 


In  Section  2.1.1,  the  velocity  potential  solution  in  cylindrical 


coordinates  is  obtained  by  evaluating  the  integral 


0°  °o  i(nxx  +  nyy) 

I  2 - 2  r  —  dnx  dny 

J  t\  -  k  3 

— OO  —CO  1  ft 


x  v 

We  first  rotate  coordinates  so  that  n  =  —  n'  -  —  n'  and 

x  r  x  r  y 

Hy  =  ^  iV  +  ^  Hy  ,  where  r^  =  x^  +  y^  .  If  we  then  drop  the  prime  (') 

on  n'  and  n'  ,  Equation  (Al)  becomes: 
x  y 


’  x 

j  T”  ~2  d\  d"y  • 

-  -  n  -  kn 


This  integral  may  be  evaluated  by  first  performing  the  integration  over 
2  2  2 

n  .  Since  n  =  1  +  n  ,  the  integrand  will  have  singularities  on 

x  x  y  _ _ 

/ ~2  2 

the  real  axis  at  n  =  +  /k.  =  n  •  These  roots  are  real  since  n 

x  —  n  y  x 

and  0^  correspond  to  the  x  and  y  components  of  the  modal  wave 
number  k^  and,  as  shown  below,  will  always  be  less  than  or  equal  to 


n  =  k  cos0 
x  n 


n  =  k  sin9 
y  n 


wavefront 
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The  integral  over  Hx  is  then 

00  irn 

**  ■  |  - - - 1 . . . . . .  •  <«> 

(n  +  /k  2  -  n  2)  (n  -  A  2  -  n  2) 

x  n  y  x  n  y 

The  integral  is  evaluation  by  integration  in  the  complex  2-plane  using 
Cauchy's  residue  theorem  and  the  contour  in  Figure  34.  In  the  figure, 
the  two  roots,  z^  and  ,  have  been  given  very  small  complex 
values  e  so  as  to  move  the  roots  away  from  the  real  axis.  In  Figure  34, 
£  is  chosen  positive  to  insure  that  the  radiation  condition  holds. 


Imz 


Figure  34.  Contour  for  nx  Integration 


The  integral  on  (-»  ,  °°)  may  then  be  broken  down  as 
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2tt  i  Z  Residues 


It  may  be  shown  that 


0  and  the  integral  is  then  equal  to  the 


j 

CR 

residue  of  the  pole  z^  enclosed  by  the  integration  contour.  The 
positive  root  is  chosen  in  closing  the  contour  to  insure  propagation 
away  from  the  source.  We  have  then  for  the  integration  over  Hx  that 
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I  =  TTi  — 


(A4) 
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when  £  goes  to  zero. 

The  integral  on  t|  to  be  evaluated  is  then 
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(A5) 
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Noting  that 
becomes 


'  2  2 

k  sinG  and  A.  -  n  *  k  cosG  ,  Equation  (A5) 
n  n  y  n 
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which  has  been  identified  by  Korn  and  Korn  (1968)  as  Sommerf eld ' s 
integral.  The  integral  is  evaluated  by  letting  0  become  complex 
with  the  integration  limits  shown.  For  these  limits,  Korn  and  Korn 
show  that 
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